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ABSTRACT 

Context. This is the second paper in a series of papers showing the results of extrasolar planet population synthesis calculations using 
our extended core accretion model. In the companion paper (Paper I), we have presented in detail the methods we use. In subsequent 
papers, we shall discuss the effect of the host star's mass on the planetary population and the influence of various properties of 
protoplanetary disks. 

Aims. In this second paper, we focus on planets orbiting solar-like stars. The goal is to use the main characteristics of the actually 
observed extrasolar planet population to derive in a statistical manner constraints on the planet formation models. 
Methods. Drawing initial conditions for our models at random from probability distributions derived as closely as possible from 
observations, we synthesize a number of planetary populations. By applying an observational detection bias appropriate for radial 
velocity surveys, we identify the potentially detectable synthetic planets. The properties of these planets are compared in quantitative 
statistical tests with the properties of a carefully selected sub-population of actually observed extrasolar planets. 
Results. We use a two dimensional Kolmogorov-Smirnov test to compare the mass-distance distributions of synthetic and observed 
planets, as well as the one dimensional version of the test to compare the M sin i, the semimajor axis and the [Fe/H] distribution. We 
find that while many combinations of parameters lead to unacceptable distributions, a number of models can account to a reasonable 
degree of statistical significance for most of the properties of the observed sample. We concurrently account for many other observed 
features, e.g. the "metallicity effect". This gives us confidence that our model captures several essential features of giant planet 
formation. In addition, the fact that many parameter combination could be rejected, indicates that planet population synthesis is 
indeed a promising approach to constrain formation models. Our simulations allow us also to extract a number of properties of the 
underlying exoplanet population that are not yet directly detectable. For example, we have derived the planetary initial mass function 
(PIMF) and have been led to conclude that the planets detected so far represent only the tip of the iceberg (9%) of all the existing 
planets. The PIMF can also be used to predict how the detectable extrasolar planet population will change as the instrumental precision 
of radial velocity surveys improves from ~ 10 m/s to ~1 m/s, or even to an extreme precision of 0.1 m/s. 

Key words. Stars: planetary systems - Stars: planetary systems: formation - Stars: planetary systems: protoplanetary disks - Planets 
and satellites: formation - Solar system: formation 



1. Introduction 

In the first paper of this series (Mordasini et al. 2008 , hereafter 
Paper I), we have presented our methods to synthesize popula- 
tions of extrasolar planets. We have explained how we use our 
extended core accretion model (Alibert et al. |2005al l to gener- 
ate synthetic planetary populations by varying in a Monte Carlo 
fashion four key variables describing the initial conditions in our 
planet formation model. As shown in Paper I, we have tried in 
deriving the probability distribution of these four variables to 
stay as close as possible to actual observations. 

We have found that the large spread of initial conditions 
resulting from the variation of the characteristics of the proto- 
planetary disk (abundance of heavy elements, mass and life- 
time) and their relative probability of occurrence leads to the 
formation of a synthetic population of planets characterized by 
a large diversity. Hence, we argued that, within the core accre- 
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tion paradigm, the observed diversity of exoplanets is a natural 
consequence of the diversity of disk properties. 

In Paper I, we have also identified a number of typical phases 
planets undergo during their formation, and found that these 
phases lead to characteristic planetary formation tracks. These 
tracks determine the final position of each planet in the distance 
to star versus planetary mass diagram (a-M) and therefore can 
be used in order to interpret the corresponding observational di- 
agram. 

Unfortunately, not all model parameters can be constrained 
by observations of proto-stellar disks. To circumvent this prob- 
lem, we present in this paper an approach that consists in com- 
paring statistically the overall characteristics of our synthetic 
planets with those of a carefully selected sub-population of actu- 
ally detected exoplanets. This approach has been made possible 
by the large number of exoplanets that have been detected over 
the recent years which has allowed to go beyond the character- 
istics of individual objects and define the characteristics of the 
ensemble population. Many studies have discussed from a obser- 
vational point of view the statistical properties of the extrasolar 
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planets, analyzing various distributions and correlations in or- 
der to address the following (and many more) issues, as recently 
reviewed by Udry & Santos ( 12007I) . 

(1) Before the detection of 51 Peg b 13 years ago by Mayor 
& Queloz ( 1995 ) it was not clear if planets outside our own Solar 
System existed, although from a theoretical point of view, there 
was no reason to doubt it. Nowadays we know that roughly 5-10 
% (e.g. Marcy et al. 120051 Cumming et al. 120081 ) of solar-like 
star in the solar neighborhood harbor a giant planet within a few 
AU in distance. 

(2) Detection biases still hinder the exploration of the full 
planetary mass domain. It is however clear that the mass distri- 
bution increases towards small mass planets (e.g. Butler et al. 
2006 1 Jorissen et al. 1200 11 1, which points towards the existence 
of a large number of yet undetected low mass planets. It is also 
known that there are very few objects with masses larger than 
~ 15 Jupiter masses inside a few AU (e.g. Marcy & Butler 2000) 
defining the "brown dwarf desert". With the detection of smaller 
and smaller mass planets, new, finer structures in the mass dis- 
tribution, like a bimodal shape at very low masses (Mayor et al. 
2009 ) have recently been suggested. 

(3) The distribution of semimajor axes consists of a pile up of 
Hot Jupiters at about 0.03 AU, followed by a relative depletion 
(the "period valley") and finally an increase in frequency further 
out at about 1 AU (e.g. Udry et al. 120031 1. Outside a few AU 
the limited time duration of the surveys does not allow definitive 
statements yet. 

(4) The combination of mass and distance has shown that 
there is an absence of massive planets at small orbital distances 
(e.g. Zucker & Mazeh 2002), and a positive correlation of plan- 
etary mass and distance (e.g. Jiang et al. 12007b . Low mass, 
Neptunian planets seem to be characterized by a different dis- 
tribution than giant planets (Udry & Santos 120071 ) . 

(5) Soon after the first discoveries of extrasolar planets it was 
noticed that the detection probability of giant planets increases 
with stellar metallicity (Gonzalez I1997I I. This "metallicity ef- 
fect" is now very well established (e.g. Fischer & Valenti 2005 1 
Santos et al. 12003b . Also correlations between stellar metallic- 
ity and the planetary semimajor axis have been discussed (e.g. 
Sozzetti 2004), but no definitive conclusions can be drawn at 
this time. The stellar mass certainly also plays a role for planet 
formation. Observations of stellar types other than FGK underly 
certain complications, but a positive correlation between stellar 
mass and frequency of massive planets seems now to be clear 
(e.g. Lovis & Mavor l2007l l. 

(6) It was found that planets in relatively tight binary sys- 
tems have statistically different properties (e.g. Eggenberger et 
al. 12004; Desidera & Barbieri 2007). For example, stars in bi- 
naries have close-in very massive planets, absent around single 
stars (e.g. Zucker & Mazeh 2002). This points toward a possible 
role of the environment in planet formation. 

(7) The observed population is now known to have an eccen- 
tricity distribution that is similar to the one of stellar binaries, 
although a group of long period, low eccentricity giant extraso- 
lar planets more similar to our giants exists (e.g. Halbwachs et 
al. l2005l 

In this work, we are addressing the points 1 to 5 while the 
remaining ones are for the moment beyond the capabilities of 
our model. In particular we wanted to investigate if by varying 
the otherwise unconstrained model parameters over a reasonable 
range of values, it was possible to reproduce, in a statistical sig- 
nificant manner, as many as possible of the characteristics of the 
observed exoplanet population listed above. To do so requires 
the ability to "observe " the synthetic planets with the same de- 



tection biases as actual observations in order to extract the di- 
rectly comparable sub-population. Since we use as the obser- 
vational comparison sample planets that have been detected by 
radial velocity (RV) methods, we have used the detection biases 
relevant for this type of planet searches with different intrinsic 
instrumental precisions. A similar approach can be used for any 
type of detection biases relevant to different detection techniques 
(transits, lensing, astrometry, direct imaging). Since these other 
techniques are sensitive to other planet characteristics than RV 
techniques, such comparisons would provide additional and in- 
dependent checks of the model which we intend to carry out 



in future work. A first such example is shown in 5 5.9.4 where 
we compare the amount of heavy elements of close-in synthetic 
planets with the one derived from internal structure modeling of 
transiting Hot Jupiters (Guillot 12008]) . 

Another distinct advantage of population synthesis calcula- 
tions is that they allow to study the global consequences of a cer- 
tain physical mechanism as shown by Ida & Lin (2008b ) for the 
example of "dead zones". In this paper, as we concentrate on the 
comparison with the observed population, we only discuss two 
such effects, namely the absence of solids inside a given semi- 
major axis and type I migration, and postpone further studies to 
subsequent publications. 

The outline of the paper is as follows: In section §|2]we sum- 
marize very shortly the methods described in Paper I we use to 
obtain the planetary populations. Section ^describes the proce- 
dure to sort out the detectable synthetic planets from the whole 
population using a synthetic observational bias. The following 
section ^4] shows how we have statistically compared this sub- 
population with the real exoplanets. The results concerning this 
comparison are given in section §|5] while those concerning the 
predictions for extremely precise radial velocity survey are in 
§(6] The conclusions are drawn in the last section, ^7] 

2. Methods 

In Paper I, we have described in detail the six step method we 
use to synthesize extrasolar planet populations. We have in par- 
ticular described the (small) changes we made to our extended 
core accretion model that was presented in Alibert et al. (2005a), 
necessary to allow for the very large number of calculations oc- 
curring for population synthesis. In our model we solve as in 
classical core accretion models (e.g. Pollack et al.[l996) the in- 
ternal structure equations for the forming giant planet, but at the 
same time we include disk evolution (using the a formalism, see 
Papaloizou & Terquem 1999) and type I and II planetary migra- 
tion. 

We have then described the four Monte Carlo variables that 
describe the varying initial conditions for planet formation, ex- 
plaining in particular how we have derived their probability dis- 
tributions from observations of (mainly) circumstellar disks. The 
random variables are (1) the dust-to-gas ratio /b/G which is con- 
strained by observed stellar metallicities [Fe/H] (Murray et al. 
1200 It Santos et al. 2003), (2) the initial gas surface density So 
at ao = 5.2 AU which is constrained by observed disk masses 
(Beckwith & Sargent 1996), (3) the rate at which photoevap- 
oration occurs M w , which determines together with a the disk 
lifetime, and is therefore constrained by observations of Haisch 
et al. (120011 ). and (4) the starting position of the embryo a stmt . 
Planetary seeds are allowed to start only in these parts of the 
disk where the isolation mass M; so is larger than the initial em- 
bryo mass M em b,o = 0.6M e and where the starting time f start of 
an embryo which is the time needed to build up such an objects 
is shorter than the disk lifetime. 
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Each population is also characterized by a number of param- 
eters that are kept constant for all planets. The most important 
(and worst constrained) parameters of the model are the viscos- 
ity parameter a for the gas disk and the efficiency factor for type 
I migration fi. By synthesizing populations using various com- 
binations of a and fi ( §5.9[ ), and comparing them to the observed 
population, we have found that the model with a — 7 x 10 3 and 
fi = 0.001 provides the overall best statistical results. We define 
this model as the nominal model. This population was already 
presented and discussed in Paper I. 

For each planet we record the full time evolution from the 
seed embryo to its final mass and position (cf. planetary for- 
mation tracks in Paper I). For the statistical analysis presented 
here, we only use the final characteristics of the planets: The fi- 
nal semimajor axis a, the total mass of accreted planetesimals 
Mhewy, the mass of the envelope M env , the total mass M, and the 
formation time of the planet. Note however that we also posses 
the fractions of icy and rocky material that were accreted dur- 
ing the formation, which allows us to study certain objects as 
GJ 436 b in detail (Figueira et al. 2008 ). To be able to compare 
quantitatively with radial velocity observations we also compute 
projected masses Msin/ for which we assume a random orien- 
tation of extrasolar planetary orbits relative to the Earth. 

3. Detection biases 

To statistically compare the characteristics of the synthetic pop- 
ulation to the actually observed one, we need to go one step 
further. We must identify which subgroup of synthetic planets 
actually could be detected in an observational survey. For this 
purpose, it is necessary to understand and quantify the various 
detection biases entering into the observational process and ap- 
ply them to the synthetic set. Such biases are technique, instru- 
ment and very likely also observer dependent. In this work, we 
have considered only the biases affecting the radial velocity (RV) 
technique, as so far the large majority of exoplanets has been de- 
tected using this technique. However, our approach can be ap- 
plied to any observational technique for which a detection prob- 
ability can be calculated as function of semimajor axis and mass 
or planetary radius. 



3.1. Synthetic RV bias 

To first order, the planet detection probability based on RV mea- 
surements increases with increasing planetary mass and decreas- 
ing distance. The instrumental precision €rv then determines 
whether the planet can be detected or not. But in fact, a large 
number of other quantities also affect the detection probability: 
The magnitude of the star, its rotation rate, the orbital eccentric- 
ity of the planet, the actual measurement schedule, stellar jitter 
and more. Using the method originally developed by Naef et 
al. (|2004 2005) for the spectrograph ELODIE, we determine 
by a x 1 analysis in a two dimensional grid in planetary period 
(1 < P < 40 000 days) and mass (1 < M < 12 720 M e ) on each, 
out of a total of 5612 grid points, the fraction of 50000 ran- 
domly chosen planetary orbits that can actually be detected by a 
spectrograph of a given precision €rv, taking into account all the 
effects mentioned above (see Naef et al. |in prep.| for a detailed 
description). This fraction represents the detection probability 
corresponding to a given planetary period and mass. 

The results of these calculations for €rv-IQ m/s are illus- 
trated in fig. [T] The graph shows the detection probability as a 
function of semimajor axis for a choice of five planetary masses 
between 100 M @ and 3000 M m . We see for example that planets 




Fig. 1. The detection probability as a function of semimajor axis 
a for five different planetary masses between 100 and 3000 M m 
and an instrumental accuracy of €rv-\Q m/s, as used for our 
synthetic RV survey. The stellar mass is 1 M Q , and it is assumed 
that the observations completely cover at least one orbital period. 



with a mass higher than about five times the mass of Jupiter Mq^_ 
(~ 1500 M ffl ) can be detected with very high probability out to 
a distance of about five AU. For a planet of about one Mq^, the 
detection probability falls below 50% outside roughly 2.5 AU. 
Note also the stroboscopic effects at orbital periods of one and 
two years, leading to a reduction of the detection probability at 
the corresponding semimajor axes. 

To decide if a synthetic planet would have been detected, 
we first determine its detection probability by interpolation in 
the detection probability grid, and then draw a random number 
between and 1 . If the latter is smaller the detection probability, 
the planet counts as one of the A^ b S synt detectable synthetic planet 
and will be used for the statistical tests (Q. 

3.2. Synthetic RV survey 

With the synthetic RV detection bias at hand, we can construct a 
synthetic RV survey by "observing" the planet population com- 
ing out of the model. To do so, we have to specify two quantities 
that characterize our MC survey: First its instrumental precision 
erv.mc, and second its temporal duration trv,mc- As we want to 
statistically compare the subgroup of detectable synthetic plan- 
ets with the real observations, these two quantities should rep- 
resent a time and instrument average over the real RV surveys 
conducted by various teams using various instruments over the 
last several years. Defining such averages is not trivial, as the in- 
strumental accuracy has changed from typically 10 m/s at the 
times of the discovery of 51 Peg b (Mayor & Queloz 1995) 
down to less than 1 m/s with HARPS (Pepe et al. 2004). This 
would have the effect that at smaller semimajor axes, planets 
of a smaller mass are known, so that 6rv,mc also should be a 
function of time, and weighted by the contributions of the dif- 
ferent observer teams. For simplicity we have assumed for our 
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synthetic MC survey a constant 6rv,mc=10 m/s and a survey du- 
ration of trv,mc=10 years. The later is needed since generally 
planet discoveries are only announced when one full orbital pe- 
riod has been covered with observations (Cumming et al. 2008 ). 
Therefore, our subgroup of potentially detectable synthetic plan- 
ets contains only planets with an orbital period smaller than 
trvmc =10 yrs corresponding to a semimajor axis of about 4.6 
AU'. 

In section ^6]we have studied the effects of changing 6rv,mc 
and trv,mc, for which we have recalculated our bias tables for 
£rv,mc=1 m / s an d 0.1 m/s. 

3.3. Observational comparison sample 

Of the more than 300 currently known extrasolar planets, not 
all can be used for quantitative comparisons with our model, 
as some planets conflict with the fundamental assumptions on 
which the model is based (Alibert et al. I2005al Paper I). We 
therefore have to isolate the ones appropriate for the statistical 
test. In particular, we select extrasolar planets based on the fol- 
lowing criteria: 

3.3.1. Sample selection criteria 

(1) The evolution of a planet may be significantly altered by the 
presence of another (massive) planet (Alibert et al. 2005b 2006; 
Thommes et al. 120081 1. In contrast, we follow the evolution of 
just one embryo per disk (Paper I), and therefore do not use any 
exoplanet that is member of a known extrasolar planetary sys- 
tem. It is clear that many single extrasolar planets could in fact 
be members of multiple systems with the small planets not de- 
tectable today. In this case, we argue that these small planets did 
not affect significantly the evolution of the massive one known 
today. 

(2) As explained in Paper I, our model doesn't describe ad- 
equately planets that migrate closer to the star than a to uch ~ 0. 1 
AU. For qualitative comparisons, we have thus to exclude all ob- 
served planets that fall into this a — M domain. This especially 
means that Hot Jupiters are excluded, reducing significantly the 
number of comparison planets. The rate of occurrence of Hot 
Jupiters is however a too important constraint on migration to 
ignore it, so we still take it into account separately, as described 
in ^6] 

(3) Planets in binary or multiple stellar systems have dif- 
ferent statistical properties than planets orbiting single stars 
(Eggenberger et al. 2004; Desidera & Barbieri 2007 ). This could 
be due to different migration mechanism (Eggenberger et al. 
2004). We therefore do not consider planets in stellar systems 
with a binary separation less than 300 AU. Binaries with a 
wider separation than ~ 300 AU have planets which exhibit no 
significant statistical differences from those around single stars 
(Desidera & Barbieri 120071 ). 

(4) In all simulations presented here, the stellar mass M* is 
fixed to 1 M Q . For stars with masses not too different from the 
solar mass which are the primary targets of RV surveys, say for 
FGK stars which have masses between 0.7 < M*/M e < 1.3, no 
evident correlation between M t and planetary properties have 
been found up to date (Udry & Santos 2007). For example, 
Fischer & Valenti (2005), hereafter FV05 don't see a correla- 
tion between stellar mass and formation probability of gas gi- 
ant planets for this small domain of primary masses. However, 
for stars with significantly different masses (M dwarfs, namely), 
theoretical studies (Laughlin et al. 120041 Ida & Lin 120051 ) in- 



dicate planetary populations with significantly different proper- 
ties, containing for example less giant planets. Indeed, observa- 
tions also indicate that giant planets are less frequent around M 
dwarfs. Ongoing dedicated M dwarf surveys e.g. with HARPS 
(Bonfils et al. 2005 ) will help to constrain this issue further from 
an observational point of view. On the other hand, intermedi- 
ate mass stars (1.5 < M*/M Q < 4) seem to have more and more 
massive giant planets (Lovis & Mayor 2007 ). Due to this reason- 
ing, we only include planets of stars with 0.7 < M*/M G < 1.3 in 
our comparison. 

(5) Many extrasolar planets have high, sometimes even very 
high eccentricities. In our model in contrast, synthetic planets 
can only have circular orbits (Paper I). Even though planet-disk 
interactions can pump eccentricity under certain circumstances 
too (Goldreich & Sari 2003 ), the high eccentricities have been 
interpreted as mainly the result of gravitational interactions be- 
tween (proto-)planets in initially more crowded systems (e.g. 
Rasio & Ford 1996 ; Adams & Laughlin 2003 ; Veras & Armitage 
2006 ). Thus, even if there is only one highly eccentric planet left 
today orbiting a certain star, its high eccentricity could be an in- 
dication that the number of planets might have been larger during 
the formation epoch, and that planet-planet interaction have been 
important in the system. Hence, planets with a high eccentricity 
are more likely than low eccentricity planets to have conflicted 
during their formation with our criterion (1). We have therefore 
excluded planets with an eccentricity larger than e max = 0.3 

(6) In the last few years, the accuracy of RV measurements 
has improved significantly, bringing the detection limit down to 
less than 1 m/s in some cases (Mayor et al. 2009). Our syn- 
thetic RV survey has in contrast a precision crv.mc = 10 m/s. We 
therefore use the same synthetic bias as for the synthetic planets 
(jj |3.1| l to sort out also these known exoplanets that could have 
been detected by our synthetic survey i.e. we only consider real 
exoplanets with a period of less than 10 years and a sufficiently 
large mass to be detectable by our synthetic bias. 

Additionally, we only include planets detected by the ra- 
dial velocity method, as other methods have a different detection 
bias. For the comparison sample of actual extrasolar planets, we 
use the compilation of observational data available online at J. 
Schneider's Extrasolar Planet Encyclopedia^] 

After applying the six criteria mentioned above, there are 
only A^bsreai = 32 known extrasolar planets left. A larger number 
would obviously be very desirable, which shows the importance 
of persistent unbiased observational campaigns. 

3.3.2. Representativity of the sample 

The low number also makes one wonder whether or not this 
small sub-sample represents well the overall sample of all known 
exoplanets. Our selection criteria in particular exclude planets 
which probably underwent strong planet-planet interactions (cri- 
teria 1 and 5). Thommes et al. (2008) have shown that disks 
which lead to the formation of just one giant gaseous planet 
might only represent a very specific class of all planet-forming 
disks, namely those near the threshold for giant planet forma- 
tion. When the mass of the disk is in contrast e.g. well above 
the threshold, prolific giant planet formation could occur. This 
could in turn be reflected in the final properties of the planets. It 
is therefore also from this perspective interesting to compare our 
32 planet sample with a larger set of exoplanets. 

To do so, we disregarded criterion 1 or 5, or both together 
and studied if the associated a - M, M, a and [Fe/H] distribu- 

1 http://exoplanet.eu 
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tions are significantly different from the iVobsreai = 32 case. We 
find that this is in general not the case. We also find that for the 
differences that do occur, the eccentricity criterion is more im- 
portant than the multiplicity criterion alone. The only difference 
that is probably significant occurs for the mass distribution, with 
the more eccentric planets being more massive. Such a corre- 
lation among the known extrasolar planets was already pointed 
out some time ago by Marcy et al. (2005 ), and can indeed be ex- 
plained by planet-planet interactions: More massive disks pro- 
duce more massive planets, and in higher numbers, leading to 
stronger scattering and therefore higher eccentricities, a behavior 
that is indeed seen in the simulations of Thommes et al. (|2008 ). 
This effect is even strengthened by the fact that in planet-planet 
interactions, the less massive bodies tend to get ejected. The vi- 
sual impression that the two distributions are however not ex- 
tremely different is confirmed by a Kolmogorov-Smirnov test 
(cf. the next section below) indicating a significance of still 37 
% that the two samples are identical. Considering the semimajor 
axis distribution, our 32 planet sample seems to contain slightly 
more planets at smaller distances (a < 1 AU). This difference 
is however clearly not statistically significant. For the metallic - 
ity, no clear trends are visible either. One can see an absence 
of very high eccentricity (> 0.6) planets at low [Fe/H]< -0.1. 
But the KS significances are always higher than 85 % (as for a) 
that the samples are identical, so this correlation is with the still 
low number of planets not significant. The general lack of sig- 
nificant correlations between a, e and [Fe/H] has already been 
found elsewhere (FV05; Udry & Santos |2D07l >. 

We conclude that in general our comparison sample, despite 
being small, represents quite well the overall giant planet popu- 
lation around solar like stars, except for a likely shift to some- 
what smaller masses. This result could indicate that at least for 
giant planets, our one-embryo-per-disk approach (see Paper I) 
leads in a statistical sense not to completely different results than 
the real multi-body formation process (Thommes et al. [2008). 
For smaller mass planets, such a generalization might however 
be less well-founded. 

4. Statistical analysis 

To assess the statistical significance of our results we perform 
four Kolmogorov-Smirnov (KS) tests in which we compare the 
distributions of the most important properties of the potentially 
detectable synthetic and the real planets. The null hypothesis is 
that the synthetic potentially detectable planets and the real plan- 
ets are drawn from the same parent distribution. Small values of 
the significance level S returned by the test show that this null 
hypothesis has to be rejected with a high probability of 1 - S . In 
particular, we perform three 1 -dimensional tests (for the distri- 
butions of M sin i, a and [Fe/H]), and one 2-dimensional KS test 
in the a - M sin i plane. 

In the 1-D case, we run standard two sided KS tests (Press 
et al. 1 1992l l. comparing the sub-population of N b SS ym detectable 
synthetic planets with the observational comparison sample con- 
taining A^obsreai = 32 real extrasolar planets. In 1-D it is possible 
to directly calculate the significance level S once the KS distance 
c/ks is known (Press et al. 1 19921 ) as 



N obsreal^ obssynt 
N obsreal obs: 



synt 



X = ( y/N e + 0.12 + 0.1 l/y/N e )xd KS 

oo 

S(x) = 2^(-iy'- 1 exp(-2/x 2 ). 



(1) 
(2) 
(3) 



In 2-D, the distribution of the KS distances c/ks for the null 
hypothesis is not independent of the shape of the distributions 
that are being examined (Press et al. [1992 ), and thus the analyt- 
ical transformation of the c/ks values to the significance S is not 
accurate in all cases. Therefore, we have proceeded as described 
by Press et al. d 1 992b in order to get the significance directly. 
After having generated a large number of detectable synthetic 
planets (iVabssynt)> that serve as the synthetic comparison popula- 
tion, we generate another A^ sam p samples of synthetic detectable 
planets, each one containing the same number of synthetic plan- 
ets A^obsreai = 32 as the real observational comparison sample. 
For each of these synthetic bootstrap samples we compute their 
KS distance c/ks to the /V bssynt planets of the detectable synthetic 
sub-population. Finally we also compute the KS distance of the 
real observations, and calculate the fraction of cases of N ydmp 
where these synthetic c/ks exceed the c/ks from the real data. This 
fraction is then the significance S . 

For each bootstrap sample we also compute as a check the 
three KS distances c/ K s for the 1-D tests with this procedure, so 
that we also get the three significances for M sin i, a and [Fe/H] 
in this direct way without using the equations above that link c/ks 
and S. As expected for the 1-D case, the two methods always 
yield very similar results. 



4.1. Defining observational constraints on the results 

One of our main goals is to test if it is possible to reproduce all 
(or at least the most important) observational characteristics at 
the same time with one single synthetic population. We consid- 
ered the following six observational characteristics as constraints 
to the model: (1) A high statistical KS significance Sa-M for the 
two dimensional distribution in the a - M sin i plane ( §5. 2} , idem 
for the one dimensional distributions of (2) the mass Sm Q5.3| >, 
(3) the semimajor axis 5 a (§5.4i, and (4) the metallicity 5[F e /H] 
(^5.5|>, then (5) a Hot Jupiter fraction Fp\ which is compatible 



with observation (4 5.6 1, and finally (6) a correct reproduction of 
the "metallicity effect", i.e. the increase of the detection proba- 
bility with stellar metallicity (jj|5.7|i. 



7=1 



We also compared the overall detection probability P of our 
synthetic survey (fraction of embryos that grew to become de- 
tectable planets) with the actual values, but we should bear in 
mind that our Monte Carlo simulations yield strictly speaking 
a different result than the observations due to the one-embryo- 
per-disk simplification: We can calculate the probability that one 
specific embryo with a given a start and r start becomes detectable, 
whereas observations yield the fraction of stars for which any 
of the initially numerous embryos in the disk finally became a 
detectable (giant) planet. 

To see if we can fulfill the observational constraints, we 
have generated several populations, keeping the probability dis- 
tributions which are constrained by observations fixed, but vary- 
ing some parameters (mainly a and fi). The varied parameters 
are listed in table [4j and their influence is discussed in |5.9| 
Typically, we were confronted with the fact that changing one 
parameter had multiple effects, bringing our results closer to one 
observational constraint, while at the same time the results for 
another deteriorated. However, many combinations resulted in 
populations that were clearly not compatible with observations. 
These negative results provide, in some sense, as much useful 
information regarding planet formation models as the positive 
ones. 
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5. Comparison with observation 



Table 1. Basic results of the population synthesis. 



5.1. Statistical assessment 

In table [TJ the basic results for the nominal population i.e. the 
population with the overall best results when compared to the 
actual population are summarized. The total number of initial 
conditions that were drawn is 7Yi n i t = 70 000. However, contrary 
to what may happen in real systems, we only start a formation 
calculation if the initial conditions are such that somewhere in 
the disk a sufficiently massive body (> M em b,o = 0.6) can form 
during the disk's lifetime (see Paper I for an explanation). 

The later conditions is fulfilled in 7V synt = 50204 (« 72 %) of 
all disks. The corresponding N syrit planets constitute what we re- 
fer to as the "full population" despite the fact that it is incomplete 
at low masses (see next). In the remaining iV noca i c = 19796 disks, 
a low fu/G coincided with a low £ so that the isolation mass in 
the disk is < 0.6 M e everywhere, and/or the disk lifetimes is 
so short (high M w together with a low Sq) that the disk disap- 
pears before such an embryo can form. For such disks which are 
hostile to planet formation in general and to giant planet forma- 
tion in particular, we don't explicitly calculate the formation of 
a planet, as we have found that detectable synthetic planets can 
form only if the disk lifetime is at least ~ 0.5 Myr (§5.8 1 and 
the isolation mass is larger than about 3 Earth masses (usually 
it is of the order of 8 M e or larger). Therefore, detectable syn- 
thetic planets cannot form in the A^ noca i c disks. Thus, the fact that 
/Vnocaic disks/initial conditions are discarded has no influence on 
the statistical analysis of the detectable sub-populations which is 
our primary interest in this paper. 

Qualitatively we expect that in the Af noca i c disks planets will 
form also. Most likely, a system of very low mass planets (less 
than a few M e ) will eventually emerge. However, as their forma- 
tion is likely to occur on timescales significantly longer than the 
gas disk lifetime, these planets will not be able to accrete nebu- 
lar gas. The resulting incompleteness of our model at low masses 
should be kept in mind when considering our predictions regard- 
ing the full population, for example in the initial mass function 
(fig# 

Within the t rv ,mc = 10 yrs survey length and with the 
£rv,mc = 10 m/s instrumental precision, 6075 planets with 
o- > Qtouch are classified as detectable in our synthetic RV survey. 
7Vhot=1386 planets migrated to a t0 uch- We call such cases "Hot" 
planets (see 5 5.6i. Assuming that all /AW planets were swal- 



Feature 



Value 



lowed by their host star, the synthetic survey has thus an overall 
detection probability of P = 6075/70000 = 8.7 %. Assuming 
the other extreme case, namely that all iVhot planets became de- 
tectable (regardless of their actual mass when they reach the 
"feeding limit" at inner boarder of our computational disk at 
fltouch ~ 0.1 AU), leads to a P whot = (6075 + 1386)/70 000 = 10.7 
%. These two values could be seen as bracketing the value 
for the real detection probability unless the effect of the one- 
embryo-per-disk approach changes the picture too dramatically 
(cf.fU). 

Nevertheless, we notice that the overall detection probabil- 
ity from the synthetic survey agrees surprisingly well with often 
quoted actual yields of 5-10 % as for example the 6.6% given 
by Marcy et al. (2005) for giant planets with a <5 AU as in our 
simulation. Cumming et al. (2008) find an extrapolated occur- 
rence rate of planets with 0.3 < M sin i/M^ < 15 of 8.5 ± 1 .3 or 
1 1 ± 1.7% out to a semimajor axis of 3 and 5 AU, respectively. 

For the KS tests, the number of bootstrap samples N samp is 
180. The specific KS results for a - M sin i, M sin i, a and [Fe/H] 
are discussed below. 



Duration of synthetic survey (t rv _mc) [yr] 10 

RV-precision of synthetic survey (6rv.mc) [m/s] 10 

Nb. of initial conditions (N m ) 70 000 

Nb. of calculated synthetic planets (N m t) 50 204 

Nb. of initial conditions without calculations (AWaic) 19 796 

Nb. of detectable synthetic planets (A^bssynt) 6075 

Nb. of synth. planets migrating to fl t0UC h (AW) 1386 

Synth, detection probability w/o AW (P) [%] 8.7 

Synth, detection probability w. AW (Pwhot) [%] 10.7 

Nb. of actual planets in obs. comp. sample (AWreai) 32 

Nb. of KS bootstrap samples (Af samp ) 180 

Significance KS a - M sin i (5 s -m) [%] 87.7 

Significance KS M sin i (S M ) [%] 95.6 

Significance KS a (S a ) [%] 63.9 

Significance KS [Fe/H] (S [Fe/H ]) [%] 21.7 



5.2. Mass-Distance diagram 

The KS test for the two dimensional distribution in the mass- 
distance plane actually checks wether planets of the correct mass 
are found at the correct distance. It is the observational con- 
straint that we weighted highest, as the mass-distance diagram 
is of similar importance for planet formation and evolution as 
the Hertzsprung-Russell diagram for stars (Ida & Lin 2004a). 
The reason for this is that it contains a lot of information about 
the planetary formation process, as described in Paper I. 

When comparing our synthetic results with the observed dis- 
tribution, one should bear in mind that our model shows the 
mass-distance distribution at the time when the gaseous disk 
disappears. Later on, it can be modified by evolutionary effects 
as evaporation or N-body interactions in initially more crowded 
systems. For the statistical comparison, we have tried to mini- 
mize those effects by carefully choosing the observational com- 
parison sample, as described in 5 3.3 



5.2.1. Full population 

In fig.|2j panel (A), the projected mass versus distance diagram 
of the full synthetic population is plotted. In the companion 
Paper I, we have thoroughly discussed this figure, so that we 
only summarize these findings here. 

One first notes that the variation of initial conditions within 
the observed boundaries results in a synthetic planet population 
with a large diversity. We conclude that the observed diversity 
of extrasolar planets is a natural consequence of the diversity of 
protoplanetary disks. Inside the envelope covered by the planets, 
various sub-structures like concentration, clumps, bars or deple- 
tions can be identified. 

The most prominent concentration can be seen near the lower 
mass boundary. There is a vast sub-population of core dominated 
low mass planets (M « M lso < 5 - 10 M e ). We call this group 
the "failed cores " in the sense that these planets did not man- 
age to grow large enough within the lifetime of the gaseous disk 
to accrete significant amount of gas. These "failed core" plan- 
ets are however not identical to the final population of terrestrial 
planets. They rather represent an earlier state in the formation 
of terrestrial planets and more massive icy counterparts beyond 
the iceline. That is to say they show the state of such planets 
at the moment when the gas disk disappears. Growth beyond 
isolation where terrestrial planets obtain their final mass occurs 
through a series of giant impacts among protoplanets of com- 
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Fig. 2. Panel (A): Projected mass Msin/ versus distance a of the full synthetic population. The feeding limit at a t ouch is plotted as 
dotted line. Planets migrating into the feeding limit have been put to 0. 1 AU. Various sub-populations which have been discussed 
in Paper I can be distinguished, namely the "failed cores", the "horizontal branch", the "main clump", and the "outer group". As 
expected, the effect of sin / is to blur slightly these structures. Panel (B): The remaining sub -population of A^ b SS ynt = 6075 actually 
detectable synthetic planets after applying the 6rv,mc = 10 m/s synthetic RV bias (small dots). The sharp cutoff at about 4.6 AU 
corresponds to an orbital period of 10 years, the assumed observational baseline t rv ,mc- The observational comparison sample with 
AWsreai = 32 real extrasolar planets is overplotted as big dots. 



parable size (which can be identified with the "failed cores") 
after the gaseous disk is gone in a final, post-oligarchic phase 
(Goldreich et al. 120041 1. Then, all "failed cores" in one disk (of 
which we however model only one) start to interact gravitation- 
ally, leading to a rearrangement through scattering and ejections 
and to mass growth by giant impacts (Ford & Chiang 2007). 
Such a behavior can be seen in the simulations of Thommes 
et al. ( 120081) . where planetary systems with a number (of order 
10) low mass planets are the typical simulation outcome for low 
mass disks. 

At semimajor axes between the feeding limit at roughly 0. 1 
AU and approximately 5 AU, and masses of 10 < M/M @ < 30, 
a "horizontal branch" of subcritical (M env /Mh e avy S= 1) cores 
is visible. Their seeds started usually outside the iceline. While 
some planets are found at intermediate semimajor axes at the 
moment when the disk disappears, the "horizontal branch" also 
also acts as the "conveyor belt" by which many low mass planets 
are transported close to the star (< 0.1 AU). Observed examples 
of the "horizontal branch" might be the Neptune -mass planets 
around HD 69830 (Lovis et al. 120061 

While migrating along the "horizontal branch", some plan- 
ets become supercritical for gas runaway accretion and leave the 
branch upwards towards higher masses (cf. the formation tracks 
in Paper I). This leads to a concentration of giant gaseous plan- 
ets at distances from the star of roughly 0.3 to 2 AU and masses 
of 100 < M/M @ < 1000. We call this concentration the "main 
clump". 

We find that once gas runaway is triggered, the planetary gas 
accretion rate becomes quickly limited by the gas accretion rate 



in the disk and the planetary mass grows larger than the local 
disk mass. As explained in Paper I, limiting the planetary gas 
accretion rate by the disk accretion rate has the consequence that 
we find a clearly less pronounced depletion of planets of inter- 
mediate masses than it is the case for the "planetary desert" 
found first by Ida & Lin (2004a). A certain depletion of plan- 
ets with masses between ~ 30 - 100 M m is however visible 
in the mass-orbit diagram of our synthetic population also as 
well as in the planetary IMF, fig. [3] as it is a typical feature of 
core accretion. The second effect (a higher planetary mass than 
the local disk gas mass), has as consequence that the migration 
mode changes from the disk dominated regime into the slower 
planet dominated type II mode (Paper I; Edgar l2Q0Tb . This slow- 
ing down prevents the planets from migrating too quickly to- 
wards the star. It also naturally causes an absence of very mas- 
sive planets at small distances from the star, as observed (Zucker 
& Mazeh|2002]). 

Planetary seeds with a large starting position a start (between 
4-7 and 20 AU), and a disk environment with a high solid sur- 
face density (high /d/g and/or So) can grow supercritical for gas 
runaway accretion in-situ and become giant planets, without the 
need of passing first though the "horizontal branch" to collect 
enough solids. This leads to the formation of another concentra- 
tion of giant planets, the "outer group " with 2 < a < 5 AU and 
1 < M/Mq^ < 20. Some of the planets in this group are thus 
so massive that they fall in the interesting category of deuterium 
burning planets (Baraffe et al. I20081 >. 

The mass-orbit diagram also shows a second depletion of 
planets at semimajor axes inside 2 AU between the "failed core" 
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planets and those in the "horizontal branch", i.e. at masses be- 
tween ~ 3 and 10 M m . At small semimajor axes (< 0.3 AU) it 
is particularly clear. The reason is the following: "Failed cores" 
which only migrate in the strongly reduced type I mode in the 
nominal case, grow only up to a mass approximately equal the 
isolation mass M iso . For the most metal rich disks considered 
here (2 = 1000 g/cm 2 , / D/G = 0.13), M iso at 0.1 AU is about 3 
M @ . Planets in the "horizontal branch" have in contrast accreted 
most of the solids in the inner part of the disk once the reach 
small semimajor axes, so that they have a minimal mass of the 
order of 10 M m at 0.1 AU. Note that the faster type I migra- 



tion, the lower this mass limit (Paper I, 5 5.9.2 1. Growth beyond 
the isolation mass up to final masses by giant impacts between 
different "failed cores" would tend to fill the depleted region. 
We can roughly estimate the mass to which Super Earth planets 
could grow by this process in-situ. If all originally present solids 
are incorporated into the planets, and their final relative spacing 
is of the order of Aa ~ a/3 as in the solar system (Goldreich et 
al. 120041) . then planets as massive as 10 Earth masses could form 
at 0. 1 AU in the most metal rich disk which would fill up at least 
partially the depleted region. 

This second depletion, which is also in visible in the plan- 
etary mass spectrum (45.3.1 1 is therefore not a robust predic- 



tion of the model, and could in principle disappear once planet 
growth after disk dispersion is included in the model. Note that 
for the statistical comparison of the detectable synthetic plan- 
ets with our comparison sample of actual known exoplanets, this 
does not constitute an issue. 



5.2.2. Detectable sub-population 

Panel (B) in fig. [2] shows the sub-population of the potentially 
detectable synthetic planets which remains after applying the 

The sharp cutoff at 



crv.mc = 10 m/s synthetic bias of §3.1 



about 4.6 AU corresponds to a 10 year period, the assumed ob- 
servational baseline trv,mc- The observational comparison sam- 
ple with A^bsreai = 32 real extrasolar planets is overplotted as big 
dots. The most striking feature is that our synthetic MC survey 
is able to detect just a small fraction of the underlying full planet 
population, between 8.7% to 10.7% of all synthetic planets (tab. 
[TJ). As expected, the planets detectable at 10 m/s are Saturn to 
Super Jupiter class planets, plus a handful planets with interme- 
diate masses (~ 50 M e ) close to the star. Even if radial velocities 
measurement in the last few years have reached a precision much 
better than 10 m/s (see § |6.2| for the detectable sub-population at 
fRV.MC = 1 or 0. 1 m/s), we still can conclude that the currently 
known extrasolar planets are just the tip of the iceberg of the real 
underlying population, as discoveries at the 1 m/s level still re- 
quire a large investment of observational time and are restricted 
to small semimajor axes. At the high mass end, we see that the 
"outer group" represents a significant reservoir of very massive 
planets at larger semimajor axes. We note that Cumming et al. 
(2008, their fig. 5) have shown that in the Keck Planet search 
program a group of very massive candidates (Msin/ > 20 M^) 
at periods > 2000 d (a > 3 AU) exists which have not yet been 
announced. Such very massive objects are virtually absent in the 
model at smaller semimajor axes (and especially do not reach 
the feeding limit), in agreement with observations (Udry et al. 
120031 . 

In the statistical comparison of the detectable sub-population 
with the observational comparison sample we find with the two 
dimensional KS test of the a - M sin i distribution a significance 
of 87.7% that the two populations come from the same parent 
population. To our knowledge, this is the first time that it is 



shown that a theoretical formation model is able to reproduce 
in a quantitative way the observed mass-distance distribution of 
an adequate comparison sample of extrasolar giant planets. 

Even if we have determined the two most important parame- 
ters of the model, a and fi by fitting the detected planet popula- 
tion, getting an agreement for any combination of parameters is 
a nontrivial result. First, a certain number of elements were given 
such as the formation model itself and the probability distribu- 
tions for fo/o, ^o, and M w which were derived from observa- 
tions. Second, the number of observational constraints that must 
be reproduced concurrently with one population is large ( §4.1 | i, 
while the number of free parameters is small. Third, at least one 
of the parameters, a, can only be varied within about one or- 
der of magnitude as observational constraints exist (King et al. 
120071 . Finally, varying parameters has complex consequences 
on the characteristics of the population thereby limiting the pos- 
sibility to force the system in a particular direction. We therefore 
interpret this result together with the other ones of this section 
as an indication that the core accretion mechanism as described 
here, while still being extremely rudimentary, must successfully 
catch several essential aspects of giant planet formation. 

5.3. Mass Msin/ 

The second distribution we have compared statistically is the 
mass distribution. It is clear that good results in the 2D a-M sin i 
distribution imply to some extent good results for the ID distri- 
butions of M sin i and a separately (whereas the opposite is not 
true). It is nevertheless worth studying these important distri- 
butions also separately, as they have been discussed extensively 
from both an observational and theoretical point of view (e.g. 
Udry & Santos 2007 ; Ida & Lin 2004b ), and because it is simpler 
to gain in this way insights into the differences between model 
and observation than in the 2D case. 



5.3.1. Planetary IMF 

Before comparing the detectable sub-population with the real 
observations, it is interesting to have a look at the underlying, 
unbiased mass distribution of the full synthetic population, as 
this can have important implications for future planet search 
campaigns. Figure [3] shows the predicted planetary initial mass 
function (PIMF) of synthetic planets around 1 M stars in the 
one-embryo-per-disk approximation. One should keep in mind 
that this PIMF shows the planetary mass spectrum at the mo- 
ment when the disk disappears. Subsequent modifications due to 
evaporation, planetary merging and ejection by N-body interac- 
tions and especially the formation of low mass planets on long 
timescales are not included in our model. 

At the largest masses, the PIMF shows a smooth decrease 
with increasing mass through the Super Jupiter (M > 3 Mo_) 
and the Deuterium Burning Planets (M > 13.6 Mq^_, cf. below in 
this section) mass domain. While this tail contains only a small 
percentage of the full population (tab. it is obviously of pri- 
mary interests for several detection techniques such as astrome- 
try or direct imaging. Planets in this domain are mainly formed 
in disks with a lifetime above average (high Eo and low M w ), 
in which the planet has plenty of time to accrete all nebular gas 
that is not photoevaporated and that viscously flows towards its 
position. 

At lower masses, in the Jovian mass domain (100 < M/M m < 
1000), the mass distribution is approximatively flat, although a 
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Fig. 3. Planetary initial mass function, corresponding to the mo- 
ment in time when the gaseous protoplanetary disk disappears. 
Several mechanisms can subsequently modify the distribution. 
The largest changes are expected to occur at low masses (below 
~ 10-20 M e ). The planetary IMF (or PIMF) has a complex 
structure with several minima and maxima (see text). 



shallow local maximum near 1-2 Jupiter mass can be seen, as it 
is observed (e.g. Jiang et al. 120071 . 

Decreasing further in mass to intermediate sized planets with 
no equivalence in the Solar System (30 < M/M @ < 100), we see 
a smooth decrease of the PIMF. It falls to a local minimum at 
about 30-40 M m . This minimum is a consequence of the same 
effect that causes the "planetary desert" (Ida & Lin 2004aJ, i.e. 
that fact that once gas runaway starts, it quickly leads to a signif- 
icant increase in mass, so that the probability for a planet to have 
a final mass just slightly larger than the one needed for runaway 
accretion is small because it is unlikely that disk dispersion cuts 
the gas supply exactly at this moment. This minimum can there- 
fore be seen as the transition between low mass, solid dominated 
planets, and gas dominated giant planets. 

It is interesting to note that this effect, which is very charac- 
teristic for the core accretion mechanism is well visible as mini- 
mum in the PIMF. Its location at a total planetary mass of about 
30 Earth masses is expected, as this is the typical total mass when 
gas runaway accretion sets in near the crossover mass, i.e. when 
M core ~ M env ~ 15 M ffi (e.g. Pollack et al. I1996I I. An obser- 
vational confirmation of the minimum at about 30 M m (which 



becomes visible at a RV precision of 1 m/s, see fig. 15 1 indeed 
seems to have occurred very recently (Mayor & Udry 2008). 
This would be a strong indication in favor of core accretion as the 
dominant giant planet formation channel, and its location would 
be a direct measurement of the mass at which runaway accretion 
typically begins. 

Note that the intermediate mass domain (30-100 M m ) is still 
quite populated in our calculations (the PIMF shows that there 
are about 2.5 times less planets in the intermediate mass than 
in the Jovian mass domain, see tab. |2j. Radial velocity surveys 
only now start to carry out observations with the precision re- 



quired to detect planets in this mass range at a significant fraction 
of an AU (Lovis et al. 12006). High precision RV measurements 
over a long time baseline will be very helpful in observationally 
characterizing this part of the PIMF from which the character- 
istic timescale for runaway gas accretion could be derived since 
the later determines how populated this intermediate mass range 
will be. This would in turn provide important information for 
formation models. For example have Miguel & Brunini (2008) 
recently shown that they find a minimum in the planetary mass 
spectrum at 100-1000 M m if they use a formation model simi- 
lar to Ida & Lin (2004a) but couple it to solid and gas accretion 
rates which fit the results of Fortier et al. (2007), rather than a 
minimum at 10-100 M s as in Ida & Lin d2004al ). 

At masses below 30 - 40 M m , the PIMF raises rapidly again 
with decreasing mass to reach a well defined local maximum 
in the Neptunian mass domain (7-30 M e ). The maximum at 
about 15 Earth masses is caused by the planets in the "horizon- 
tal branch" in the a - Msin/ diagram, i.e. of subcritical cores 
migrating in disk dominated type II and collecting solids (Paper 
I). These planets have properties similar to the ice giants of our 
own planetary system with 0.01 < M env /M core < 0.4. 

A next local minimum occurs at about 7 M e . It marks the 
boundary between "failed cores" and planets in the "horizontal 
branch". The reason for this minimum is clear as well. In or- 
der to best reproduce the observed population of giant extrasolar 
planets, especially the semimajor axis distribution, we found that 
very small type I migration efficiency factors are needed (see 
section 5.9.2 1. With such low type I migration rates, isolation 



phenomena in which the planet almost completely empties its 
feeding zone of solids at its starting position become important. 
This results in the quenching of the embryo's growth at a mass of 
roughly < 5 M @ , which approximately corresponds to the maxi- 
mal mass to which a "failed core" can grow in a disk with a mean 
solid surface density just beyond the iceline. Type II migration 
is in contrast calculated at its nominal rate. Therefore, as soon 
as the migration mode switches from type I to type II, the core 
moves quickly into new, undisturbed regions of the disk with a 
large supply of new planetesimals to accrete and therefore re- 
sumes its growth to reach masses of order M ~ 20 M m which is 
the typical mass an embryo reaches after it has swept significant 
parts of the inner regions of the disk while migrating through the 
"horizontal branch". Hence, except for unlikely timing effects, 
embryos with a final mass between the two values mentioned 
above are not as frequent as the others. 

Finally, the PIMF begins to raise very quickly in the Super 
Earth mass domain, indicative of the large population of low 
mass "failed cores" planets already seen in the mass-distance 
diagram. 

The quantitative characteristics of the planetary IMF at such 
low masses {i.e. below 10-20 Earth masses) predicted by our 
model should be considered with great caution. As we have al- 
ready mentioned, is our model essentially a giant planet forma- 
tion model and is therefore severely incomplete at the low mass 
end (see also Paper I). For example, the minimum at 7 M e be- 
tween Super Earth and Neptunian planets might be an artifact re- 
sulting from not having considered the merging of "failed core" 
after disk dispersal as well as planetary formation in the N noc . d i c 
solid-poor disks. Despite this caution, we nevertheless believe 
that a strong rise of the PIMF somewhere below ~ 10-20 M e 
is a real feature simply reflecting the fact that the vast majority 
of embryos are located in disks that do not allow the formation 
of massive cores. 

Forming Jovian planets is in contrast only possible when ini- 
tial conditions are, by chance, planet formation friendly which 
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Fig. 4. Statistical test of the mass distribution. Panel (A): Histogram of the projected mass Msini of the sub-population of the 
AWssynt = 6075 detectable synthetic planets (solid line) and the observational comparison sample of N b srea \ = 32 planets (dotted 
line). Panel (B): Cumulative distribution function corresponding to (A), showing that the two distributions are very similar. 



Table 2. Percentage of various planet types forming from the 
Mnit = 70 000 initial conditions that were generated ("Init"), of 
the detectable synthetic sub-population ("Obssynt", N b SS ym = 
6075 planets), and of the observational comparison sample 
("Obsreal", iV bsreai = 32 planets). For the latter, Poisson errors 
are assumed. For the full population, the true mass M is used, 
for the other two populations the projected mass M sin i to al- 
low comparison. "Not calculated" corresponds to the A^ noca i c = 
19796 initial conditions where no 0.6 M e seed could be formed 
during the disk lifetime. DBP stands for Deuterium Burning 
Planets. 



Type [%] 


Range [M e ] 


Init 


Obssynt 


Obsreal 


Not calculated 




28.3 






Super Earth 


< 7 


41.4 








Neptunian 


7-30 


11.8 








Intermediate 


30-100 


4.2 


2.4 


6.3+4.4 


Jovian 


100-1000 


10.4 


66.6 


68.8 ± 14.7 


Super Jupiter 


1000-4323 


3.5 


28.1 


24.9 ± 8.8 


DBP 


>4323 


0.4 


2.9 






corresponds to a disk with significant amounts of solids and gas. 
It is interesting to note in this context that the Sun is in the low 
metallicity tail of planet host stars (Udry & Santos 2007). Under 
these favorable conditions, the tip of the iceberg of the full plan- 
etary population is formed, the one we detect with radial velocity 
techniques at a precision of 10 m/s. Usually, our attention is fo- 
cussed on stars which do have detectable (giant) planets. But the 
finding that 90% of the stars in RV surveys do not appear to have 
detectable giant planets is an important result as well and is in 
full agreement with our PIMF. 



From the consideration of the PIMF, we compute the frac- 
tions of different types of synthetic planets in table [2] The num- 
bers for the lowest mass bins must very likely be seen as lower 
limits. For example, many of the Alcaic initial conditions for 
which we did not calculate the formation of a planet as no 0.6 
M @ seed could be formed during the lifetime of the gas disk, may 
eventually lead to the formation of planets with masses below 7 
M m . If this were to be the case, this mass bin would eventually 
contain 69.7 % of all initial conditions. 



It is clear that the exact value of the masses defining the var- 
ious bins in tab. [2] are somewhat arbitrary but are nevertheless 
based on the various features of the PIMF. In any case it is inter- 
esting to note that the maxima of the PIMF correspond roughly 
to planetary masses occurring in our solar system. 

Finally, one notes that the core accretion mechanism is able 
(under certain assumptions, see Paper I) to produce planets 
which lie in a mass domain where, at least in the absence of 
a solid core, deuterium fusion occurs (M > 0.012 - 0.013 
M * 12.6 - 13.6 M % , Chabrier & Baraffe 120001 . Therefore, 
we simply call all planets with a mass larger than 13.6 Mq, 
"Deuterium Burning Planets" (DBP, see Baraffe et al. [2008 ). 
Such massive planets are however rare objects: only 0.4 % of 
all disks produce a planet larger than 13.6 M^, and just 11 out 
of 70 000 initial conditions produce an object more massive than 
30 Mo_. The existence of the "brown dwarf desert" is therefore 
not in disagreement with these findings. For example, Marcy 
& Butler (2000]! have estimated a frequency of brown dwarfs 
within 3 AU of < 0.5%. The overall largest synthetic planet has 
a mass of 38 M^, which is just in the middle of the desert (Lovis 
& Mayor|2007]l. 
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5.3.2. Observed mass distribution 

In fig. [4j the two panels illustrate the results of the statistical 
comparison of the Msini distribution of the detectable sub- 
populations with the observational comparison sample. Panel 
(A) of fig. [4] shows the mass histogram of the sub-population of 
the detectable synthetic planets, together with the observational 
comparison sample. Both distributions peak at about 2 M^, The 
decrease at high masses simply reflects the decrease in the unbi- 
ased PIMF in the Super Jupiter and DBP domain. The decrease 
at low masses is mainly an effect of the observational detection 
limit of 10 m/s, but not entirely, because also the underlying dis- 
tribution decreases below ~ 1 Mq^_ . 

Panel (B) shows the corresponding cumulative distribution 
functions. The plots make it obvious that the two distributions 
are very similar. Some differences exists at the low and high 
mass end. At the high mass end this could be related to the fact 
that the -/V bsi- ea i = 32 observational comparison sample only in- 
cludes companions if their mass is smaller than 20 Mo_, while 
we do not use such a criterion, so that we might should add some 
brown dwarf candidates, as suggested by Lovis & Mayor (2007). 

Applying the KS test results in a high significance that both 
the observed and the synthetic population are drawn from the 
same parent distribution. The bootstrap method leads to 5 m = 
95.6%, while using eqs. [T]to[3]leads to 96.8 %. We have found 
that it is possible to obtain good, or at least non-zero statistical 
significance for the mass distribution even if some parameters 
were changed (cf. 5 5.9 1, in contrast to S a _M or S a , that often fall 



to virtually zero. These surprisingly good results for the mass 
distribution has led us to conclude that gas accretion by a core 
may be relatively well understood as compared to, for example, 
the migration of it. 

The good agreement with observational data is also visi- 
ble in fig. [5] where the mass distribution of the detectable sub- 
population is compared in linear bins to Marcy's et al. (2005 ) 
power law dN/dM oc M~' 05 inferred from the Keck, Lick & 
AAT observational data (See Udry & Santos 120071 for an up- 
dated version of this plot.). Similarly to the observational result 
this power law gives a good representation of the mass distribu- 
tion in the Jovian planet regime. This is due to the approxima- 
tively fiat part (in logarithmic units) of the PIMF in the Jovian 
mass domain discussed above. At masses much larger than one 
Jupiter mass (> 7 M^), and especially in the DBP domain, 
the number of planets is clearly lower than inferred from the 
dN/dM oc M~ L05 law, again in good agreement with Marcy et 
al. ((2005 ) or Butler et al. (12006b . This is due to the decrease of 
the PIMF in the Super Jupiter and DBP tail. From the complex 
shape of the unbiased PIMF in fig.[3]it is however also clear that 
one power law cannot be used to describe the full domain of 
planetary masses as already noted by Mayor & Udry (2008). 

5.4. Semi-major axis a 

The third distribution we have compared with the observations 
is the semimajor axis. In the comparison, the "Hot" synthetic 
planets that reach the inner boundary of the computational disk 
at atouch are again not included, as the structure of the disk closer 
than ~ 0.1 AU to the star is likely to be more complicated than 
described in the model, and probably causes particular effects, 
like stopping planetary migration. Correspondingly, observed 
Hot Jupiters are neither included ( §3.3). T he constraints these 
planets put on the model are treated in jj |5.6| 

In fig. |6j Panel (A) the distributions of the semimajor axes 
of the detectable synthetic planets and the A^ b sr eai = 32 obser- 
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Fig. 5. Histogram of Msini of the 6075 detectable synthetic 
planets compared to the observationally inferred dN/dM oc 
AT 1 05 power law of Marcy et al. (120051 1 . 



vational comparison sample is plotted. Although the two his- 
tograms are similar (an approximately flat distribution at smaller 
distances followed by a sharp upturn at larger distances), two 
differences exist. First, among the synthetic population there are 
more detectable planets inside ~ 1 AU than in the observational 
sample, and second, the upturn occurs in the model at a larger 
distance (~ 2 AU) instead of about 1 AU as in the observational 
comparison sample, or as e.g. also shown in Marcy et al. (2005). 

In the model, the inner fiat part is populated by giant planets 
of the "main clump" which first migrate through the "horizontal 
branch" to collect solids (Paper I), while the upturn at ~ 2 AU 
is caused by "outer group" planets which grow supercritical for 
gas runaway accretion in-situ. 

Giant planets at smaller distances thus have a different, more 
complicated formation and especially migration history, than 
those further away. The number of giant planets inside ~ 1 AU 
is thus dependent on the efficiency with which migrating cores 
can accrete planetesimals, or the solid surface density profile be- 
low 1 AU. Both these factors are described only roughly in the 
model as the accretion rate is the same for a migrating as for a 
non migrating planet, and we use a simple oc r~ 3 ^ 2 law for the 
solid surface density, which could be modified for example by 
planetesimal drift due to gas drag which occurs on the fastest 
timescales at small distances (e.g. Chambers 2006). The loca- 
tion of the "outer group" and thus the upturn at about 2 AU is 
dependent on the location of the iceline. In tests where we have 
arbitrarily reduced the location of the iceline as obtained by our 
a disk model (Paper I) by 1 AU lead to a inward shift of the 
"outer group" by approximately also 1 AU, bringing it actually 
to better agreement with the observational data. 

Panel (B) shows the cumulative distribution functions cor- 
responding to (A). The KS test leads to a significance of S. d = 
63.9% using the bootstrap samples, eqs. [T] to [3] to 63.5 %. This 
confirms the visual impression that the observed and the syn- 
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Fig. 6. Statistical test of the semi-major axis distribution. Panel (A): Distribution of the final semimajor axes of the detectable 
synthetic planets outside the feeding limit (solid line) and the observational comparison sample (dotted line). Panel (B): Cumulative 
distribution function corresponding to (A). 



thetic semimajor axis distributions differ more than the mass dis- 
tributions, even if they are still statistically similar. 

When comparing the semimajor axis distribution with the 
observations, one should keep in mind that this distribution is 
very likely more affected by long timescale planet-planet inter- 
actions (e.g. Ford & Chiang |2007l Veras & Armitage 2004) than 
the mass distribution because scattering and ejections are more 
likely than collisions (Ford & Chiang 2007), so that it is not nec- 
essarily expected that our model can accurately reproduce the 
observations. It would therefore be interesting to couple the out- 
come of a population synthesis calculations to such numerical 
N-body scattering experiments, similar to the work of Thommes 
et al. d2008l . 



5.5. Metallicity [Fe/H] 

The fourth distribution we have statistically compared is the 
metallicity distribution. The so called "metallicity effect", i.e. 
the increase of the detection probability with stellar metallicity 
is discussed in 35.71 



It is clear that such a test only makes sense if one can assume 
that the planets in the observational comparison sample origi- 
nate from a search within a sample with a similar metallicity 
distribution as the CORALIE sample. This is approximately the 
case for the 1040 FGK stars in the Keck, Lick and AAT planet 
search sample (FV05, Paper I). It is clearly not the case for de- 
tections coming from metallicity biased search programs. This 
is of course also the case for all other distribution we have com- 
pared statistically. As most planets in our observational compar- 
ison sample have rather long periods (and are discoveries dat- 
ing several years back), we assume that our statistical compar- 
ison sample is in that sense not "contaminated" by planets of a 



metallicity biased search program (see Santos et al. 2005 for a 
discussion of this point). 

Panel (A) of fig.[7]displays three distributions: The metallic- 
ity distribution of the synthetic detectable planets, the metallic- 
ity distribution of the observational comparison sample, and the 
fit to the metallicity distribution of the CORALIE planet search 
sample (Santos et al. |2003t . This distribution is a Gaussian with 
fi - -0.02 and cr = 0.22 and is the distribution from which we 
drew [Fe/H] (Paper I), and thus represents the metallicity distri- 
bution of all the A^i n i t initial conditions. 

It is an observationally very well established fact that planet 
host stars have a clearly higher mean metallicity than the com- 
plete search sample. The offset between the two distributions is 
observationally found to be 0.1-0.25 dex (Santos et al. [2004b; 
FV05). Panel (A) shows that our simulations reproduce this ob- 
servational constraint. Here too, the detectable sub-population is 
richer in metals than the full population. We find that the de- 
tectable sub-populations has a mean [Fe/H] of 0.09 i.e. we find 
a relative shift of 0. 1 1 dex. 

The corresponding cumulative distributions in panel (B) 
show that the distributions of the actual and synthetic detectable 
planets are similar. The -/V b sre ai = 32 observational comparison 
sample might be affected by small number effects, as indicated 
by the bumpy structure of the histogram. It also seems that our 
synthetic distribution suffers from a certain deficit of detectable 
planets in the clearly subsolar metallicity regime, while we find 
more planets at around [Fe/H]=0, where a decrease in the obser- 
vational data exists. This decrease is however very likely only a 
small number effect, as indicated by the comparison with a less 
specific set of exoplanets (5 3.3i. The KS test returns a clearly 
non-zero, but rather low significance for both samples to come 
from the same parent distribution of 21.7 % using the bootstrap 
method. Using eqs.[T]to[3]leads to 23.5%. 
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Fig. 7. Statistical test of the metallicity distribution. Panel (A): Solid line: [Fe/H] histogram of the detectable synthetic planets. 
Dotted line: Distribution in the observational comparison sample. Dashed line: Fit to the metallicity distribution in the CORALIE 
planet search sample, from which the A^ n i t initial conditions are drawn. The distributions of both the synthetic and actual detectable 
planets are shifted towards higher metallicities relative to this curve. Panel (B): Corresponding cumulative distribution functions of 
the synthetic and actual planets. 



We plan to include in future a self-consistent calculation of 
the early solid surface density evolution, similar to the model of 
Rozyczka et al. (120041 1 as these authors have shown that this is 
important for reproducing the observed metallicity distribution 
of the extrasolar planets. 

5.6. Observational constraints for the "Hot" planets 

As explained in Paper I, due to possible effects like Roche lobe 
overflow to the star (Trilling et al. I1998K evaporation or partial 
accretion of gas streaming past the planet onto the star we do not 
really know the fate of the "Hot" synthetic planets migrating to 
atouch ~ 0.1 AU. 

Up to this point, "Hot" planets were therefore excluded from 
the statistical analysis. However, Hot Jupiters represent an im- 
portant feature of the actual population of extrasolar planets, 
and therefore we still would like to use them to test our mod- 
els. To be able to do this, we proceeded in the following way: 
FV05 have published a list of 850 FGK stars that have enough 
RV observations that every planet with a radial velocity semi- 
amplitude K >30 m/s and a period T shorter than 4 yr was de- 
tected (uniform detectability criterion). Of the 850 stars, 47 have 
a planet (or several planets) detectable with this observational 
bias. Using the planetary orbital parameter from FV05, we have 
then calculated how many of the 47 stars with planets have a 
companion (using the most massive one if there is more than 
one planet) that would have been classified as a "Hot" planet 
in our simulations, using the a touc h criterion. This is the case 
for 9 stars. We can use this to define a simple new constraint 
in the following way. With a K >30 m/s, T < 4 yr pass-fail 
bias, the overall detection probability excluding the "Hot" plan- 
ets f FV is 38/850 ~ 4.5%. Including the "Hot" planets leads to 
a fW.whot = 47/850 ~ 5.5%. The fraction of "Hot" detectable 
planets of all detectable planets is Fpy = 9/47 ~ 19%. This last 



figure is clearly the most important observational constraint, as 
it is a ratio between detectable planets only, which is likely to 
reduce possible unwanted consequences of the one-embryo-per- 
disk approach. We also note that the fraction of all stars with a 
detectable "Hot" giant planet is 9/850 ~ 1% as in Marcy et al. 
(|2005]l who find 1.3 + 0.3%. 

We then proceeded in the same way for the synthetic plan- 
ets. We sorted out the subset of synthetic planets which have 
K >30 m/s and T < 4 yr ("FV05 bias"). The radial velocity 
semi-amplitude K and T are calculated as (Udry 2000) 
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where the second lines apply for our special case of a circular 
orbit and M, = 1M . 

Such a pass-fail detection criterion is much more rudimen- 
tary than the synthetic observational detection probability pre- 
sented in jj |3.1| but we want to use a detection criterion as simi- 
lar as possible to FV05. For the "Hot" planets, whose real final 
semimajor axis is unknown we once calculated K assuming that 
their final semimajor axis is 0.1 AU, and once assuming 0.01 
AU, but always using the mass when they arrive at the feeding 
limit at a t0U ch- The smaller assumed final semimajor axis is of 
course extremely close to the star (the radius of the sun is about 
0.005 AU). But in this way, we obtain an approximative upper 
and lower boundary. 
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Table 3. Constraints on the "Hot" planets from Fischer & Valenti (2005 ) and results from the synthetic population. The two columns 
0.1 and 0.01 AU represent the results obtained if such final semimajor axes for the synthetic "Hot" planets are assumed. 



Feature 


0.1 AU 0.01 AU 


obs. const. (FV05) 


Nb. of detectable synth. planets with the FV05 bias outside feeding limit 




3382 






Detection probability with FV05 bias excluding "Hot" planets (/Vv) [%] 




4.8 


4.5 ± 


0.7 


Nb. of "Hot" synth. planets detectable with the FV05 bias 


259 


407 






Detection probability with FV05 bias including "Hot" planets (/Vv.whot) [%] 


5.2 


5.4 


> 5.5 


±0.8 


Fraction of "Hot" detectable planets of all FV05 detectable planets (Fpv) [%] 


7.1 


10.7 


> 19 


+ 6 



The results are as follows (tab. 3}: Outside the feeding limit, 
3382 synthetic planets are detected with the FV05 bias. This re- 
sults in a P FV of 3382/70000 ~ 4.8% for the synthetic popula- 
tion, compatible with the observed value of 4.5 ± 0.7%, assum- 
ing for the error that the observations follow Poisson statistics 
(as FV05). Among the "Hot" planets, 259 or 407 planets are 
detected with the FV05 bias, for a = 0.1 or 0.01 AU, respec- 
tively. This gives us a synthetic P F v,whot of (259+3383)/70000 
to (407+3383)/70 000 corresponding to 5.2 to 5.4 %, again com- 
patible with the observations (5.5 + 0.8%). For the most impor- 
tant constraint, the fraction of detectable "Hot" planets among all 
planets detectable with the FV05 bias, we find for the synthetic 
population Fpy =7.1 - 10.7 % for an assumed final semimajor 
axis of 0.1 and 0.01 AU, respectively. While the value using 0.1 
AU is too low, the upper limit is not too far from the observed ra- 
tio of > 19+6% given the large error bars. We note that this result 
corresponds to an overall frequency of synthetic Hot Jupiters of 
259/70000=0.37 % and 407/70000=0.58 % again for the lower 
and the upper limit, respectively. 



0.3 



We note that in contrast to our result Ida & Lin (2004a ) have 
found that a much higher number of planets must have migrated 
close to the star compared to what is observed today which 
would correspond to a significantly higher F FV . Such a higher 
Fpv is also in agreement with the observational constraint, as 
planets may have perished by falling into the star or by evapo- 
ration. However, we have found in our calculations that we can- 
not obtain synthetic populations with a higher Fpy (as obtained 
when assuming a higher type I migration rate, 5 5.9.2| i without 
simultaneously increasing also the number of detectable planets 
at a < 1 AU but still outside the feeding limit. Since the nominal 
model tends to already overestimate the number of planets in this 
region, we conclude that it is not possible to increase F F y with- 
out degrading significantly the KS test result for the semimajor 
axis. This is an example that illustrates clearly that the model 
cannot be tweaked in a particular direction by a suitable choice 
of parameters. For the moment, we can only speculate that these 
difficulties are due to an incomplete description of migration, 
and/or that our simulations lack a mechanism that "produces" 
Hot Jupiters without changing the semimajor axis distribution 
at larger distances. Note that a simple stopping mechanism is 
not a solution since we actually have kept all planets migrat- 
ing to these short distances. A mechanism that is controlled by 
the thermodynamic structure of the disk that could have such an 
effect is presented in jj |5.9.5 Alternatively, one could also imag- 
ine that the very strong instrumental (and very likely also strong 
observer's psychological) bias to find Hot Jupiters is not com- 
pletely corrected for in the observational comparison data. 

Indeed, other surveys have found a rather lower rate of occur- 
rence of Hot Jupiters (less than 1 %), which we reproduce better 
with our simulations (0.4-0.6 %): Naef et al. d20051 > find for the 
ELODIE planet search campaign that a fraction of 0.7 ± 0.5% 
of stars have a giant planet with a period less than 5 days and a 
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Fig. 8. Fraction of embryos that become giant planets detectable 
with K > 30 m/s and T < 4 yrs as a function of metallicity 
for the synthetic and actual surveys. The solid line assumes that 
the "Hot" synthetic planets stop at 0.01 AU, the dashed one that 
they do so at 0.1 AU. The dotted lines are the observationally 
determined values and error bars taken from Fischer & Valenti 
(2005). The dashed-dotted line is the averaged result from the 
CORALIE and Lick-Keck-AAT samples (Udry & Santos l2007l 
The smooth solid curve is given as 0.05 x 10 [Fe/H1 i.e. oc Z, which 
approximately fits our results, while the smooth dotted curve is 
the 0.03 x lO^mFe/m i e K Z 2 fit from Marcy et al J20Q5] ) 



mass larger than 0.42 M^. Cumming et al. (2008) find a similar 
result of 0.65 ± 0.4%. Fressin et al. (120071 ) find that a fraction of 
1/215 = 0.47 % of late main-sequence stars are orbited by a giant 
planet with a period between 1-5 days, using combined results 
from radial- velocity and photometric transit surveys. 

5. 7. Metallicity effect 

The strong correlation between stellar metallicity, and the like- 
lihood to find giant planets, the so called "metallicity effect", 
is possibly the best established observational constraint (e.g. 
Gonzales [T9971 Santos et al. I2004bl Fischer & Valenti 2005). 
Even though the role of metallicity in planet formation will be 
discussed in detail in a forthcoming publication, we check here 
if our synthetic population presents the same effect as has been 
done for other giant planet formation models (Ida & Lin 2004b ; 
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Kornet et al. [2005 ; Robinson et al. [20061 Matsuo et al. 2007'). 



In 5 5.5 we have already compared the [Fe/H] distribution of the 
synthetic and the actual planets. But as explained there, these 
distributions depend on the metallicity distribution of the planet 
search sample. By dividing the [Fe/H] distribution of the planet 
host stars by the one of the full sample we can correct for that. 

To do so, we use the sub-population of 3641 respectively 
3789 planets detectable with the FV05 bias described in the last 
section and assigned them to the same bins in metallicity 0.1 
dex wide as FV05. We then normalize by dividing the number 
in each bin by the total number of initial conditions in the corre- 
sponding bin. With this procedure, both synthetic and observed 
data have been binned in identical manner which allows for an 
accurate comparison. 

Fig. [8] shows the results of this comparison. It can be seen 
that the synthetic population reproduces the observed positive 
correlation between metallicity and the detection probability of 
giant planets. Even quantitatively it agrees with the observational 
constraint from FV05 within the error bars in all bins except two 
at 0. 1-0.2 dex. Note that we have not additionally normalized the 
distribution at any bins. The fact that we also reproduce the ab- 
solute numbers well means that our single embryo is obviously a 
good representative for the large number of all nascent embryos. 

Fig. [8] makes however also clear that the detection proba- 
bility as a function of [Fe/H] in the model and the observational 
data are not identical. At largely supersolar metallicities, our pre- 
dicted fractions are below the observed values (but still within 
the Poisson error given by FV05) while at mildly supersolar val- 
ues we find slightly higher fractions than FV05. We find that 
the increase of the detection probability scales approximately 
as oc Z (Z being the stellar mass fraction of heavy elements), 
while observational data rather indicate a quadratic dependence 
on Z (Marcy et al. 2005 ), or two different regimes for metallici- 
ties larger respectively smaller than 0.0 dex, with no dependence 
on [Fe/H] in the subsolar regime (Santos et al. 2004b; Udry & 
Santos 2007 ). Our theoretical predictions seem to follow approx- 
imately the same dependence on [Fe/H] over the entire domain, 
which is the also the case for Ida & Lin (I2004bll2008bl . 

The reasons for this difference in slope could be metallic- 
ity effects that are not included in the model. At the moment, 
the only quantity where /d/g enters is the solid surface density. 
Metallicity effects that are not included are for example changes 
in disk and envelope opacity, or a change in relevant planetesi- 
mal size and consequent radial distribution. This last effect was 
found by Kornet et al. (120051 ) to increase the detectable fraction 
especially at low [Fe/H]. 

In any case, it is apparent that the synthetic population does 
quite well when it comes to reproducing the "metallicity effect". 
The fact that this happens concurrently with all other observa- 
tional constraints, is the critical point in this study. We thus con- 
firm the result of parameterized core accretion models (Matsuo 
et al. 2007 ) that the core accretion paradigm seems to be capable 
of explaining the vast majority of known extrasolar giant planets. 

A closer look at the population reveals that the metallicity 
effect is due to the combination of two effects: First, metal rich 
systems favor the formation of massive planets, as one would in- 
tuitively expect from the core accretion theory, as for metal rich 
disks the core reaches the critical mass for runaway gas accre- 
tion earlier, therefore allowing more gas to be accreted before 
the disappearance of the disk. The second is simply the fact that 
the RV technique discovers massive planets more easily (fig.[TJ. 
It does however not imply an absence of other types of planets 
(not Jovian) at low metallicity. 
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Fig. 9. Distribution of the disk lifetimes in which detectable syn- 
thetic planets are formed, which are upper limits for the forma- 
tion timescales of the detectable (giant) synthetic planets them- 
selves. Detectable planets originate from disks that are in the 
mean somewhat longer living than average disks, but giant plan- 
ets can also form in disks that have a lifetime as short as 1-2 
Myr. 



With the metallicity effect, we have treated all six observa- 
tional constraints mentioned in jj |4. 1[ In the next section, we ad- 
dress an observational constraint that is not posed by the extra- 
solar planet population itself, but by the protoplanetary disk they 
are formed in. 



5.8. Formation timescales 



Early core accretion models (e.g. Pollack et al. 1996) suffered 
from the so called timescale problem. The predicted formation 
timescales for giant planets were comparable to, or even longer 
than typically observed disk dispersion timescales. This lead to 
the hypothesis that a faster formation mechanism is needed to 
form Jovian planets, namely the gravitational instability model 
(Boss l200Tl Maver et al. l2U04l . 

Later extensions and improvements of the core accretion 
model, such as the inclusion of planetary migration and con- 
current disk evolution (Alibert et al.[2004 2005a) or the use of 
modified grain opacities in the envelope (Podolak 2003 1 Hubickj 
et al. 120051 ). showed that the core accretion mechanism is able 
to form giant gaseous planets well within the timescale limit 
imposed by the observations. On the other hand, Fortier et al. 
(120071 ) have recently shown that when they employ a core ac- 
cretion rate based on the calculations by Thommes et al. (2003 ) 
for the oligarchic growth regime, they find for identical initial 
conditions a formation time for Jupiter about one order of mag- 
nitude larger than Pollack et al. (1996), so that this question is 
not yet settled (see Paper I). In the population synthesis calcula- 
tions presented here which use the same (faster) accretion rates 
of solids as Pollack et al. (1996), the formation timescales of 
giant planets are by construction compatible with the observed 
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Table 4. Important parameters for the simulations. If several 
were tested, the value for the nominal case which gives the best 
results in the statistical analysis is printed in italic letters. 



Feature 


Values 


Type I migration efficiency factor jj 


0.001, 0.01, 0.1, 1 


Viscosity parameter a 


0.001, 0.005, 0.007, 0.01 


Initial gas disk profile E(a, t = 0) 


oc a- 3 ' 2 


Rockline included 


yes, no 


Iceline included 


yes 


Photoevaporation included 


yes, no 


/d/g,o 


0.016, 0.02 0.04, 0.05 


Distribution for X 


Taurus, Ophiuchus 


Host star mass 


1M 


Initial embryo mass M emb0 


0.1,0.3, 0.6 M e 



disk lifetimes, as the distribution of the photoevaporation rate 
M w was adjusted to reproduce the observed disk lifetime distri- 
bution (Paper I), and disk lifetimes are obviously upper limits for 
the formation timescales of gaseous planets (as we neglect any 
formation after the gas disk dispersion, the disk lifetimes are in 
fact an upper limit for the formation timescales of all planets in 
the model). 

A posteriori, the Monte Carlo simulations yield the distri- 
bution of lifetimes of those disks which eventually produced a 
giant planet in the detectable sub-population. Fig. [9] shows this 
distribution. We find that giant planets form in disks with a mean 
lifetime that is about one Myr longer than the mean lifetimes of 
the disks of all N synt planets (mean lifetimes of 3.50 and 2.49 
Myr, respectively). Some giant planets (~ 14.5%) were however 
also formed in short-lived disks (< 2 Myr). This provides a hint 
that from the disk lifetimes aspect alone, giant planet formation 
is not be completely inhibited (albeit less likely) in a dense stel- 
lar environment in which circumstellar disk dispersion is rapid. 
However, for this to occur, high amounts of solids are required: 
The mean metallicity of the disks in which giant planets form 
within 2 Myr is 0.18 dex, twice as large as the mean metallicity 
of all detectable giant planets (0.09 dex, as mentioned in 5 5.5 1. 



5.9. Influence of parameters 

The results presented up to this point relate to the nominal popu- 
lation. This nominal population was defined as the population 
that provided an overall best match to all the tests described 
above from a large set of populations that were obtained from 
varying parameters. Table [4] summarizes the values of some pa- 
rameters and distributions that were used to explore all space 
available. The synthetic populations obtained with these differ- 
ent parameters were compared to the actual exoplanets using the 
same statistical methods as discussed above. In this section, we 
briefly discuss some outcomes of varying various free parame- 
ters. Our main focus rests on determining the influence of chang- 
ing the type I migration efficiency factor fi since this has the 
largest influence on the sub-population of "Hot" planets. 



5.9.1. Varying S and M embj0 

In all cases mentioned here, only one parameter was changed 
from its default value at one time. 

Using the disk mass distribution derived from observa- 
tional data of the Taurus-Auriga star forming region (with = 
-1.66, cr = 0.74, cf. Paper I) instead of Ophiuchus for So does 
not have a marked influence on our results. In the KS tests, val- 



ues of 65, 85, 58 and 50% are found for S a -M, Sm, S a an d S [Fe/H], 
respectively. The reason for these rather small changes is par- 
tially linked to the fact that we kept the range of possible Eo 
values constant (50-1000 g/cm 2 ), so that in both cases a rather 
limited part of the full distribution is covered (Paper I). 

We have also synthesized a population with M em b,o = 0.3 
instead of 0.6 M e . As the starting mass has been chosen arbitrar- 
ily, changing it within certain limits should not have a significant 
influence on the final results. A different M &m \,fl changes the dis- 
tributions of a start (via the isolation mass criterion) and naturally 
also f start . In addition lowering M em b.o increases the number of 
disks that can form such embryos within their lifetime. The frac- 
tion of initial conditions A^ noC aic/Mnit drops from 28% to 21%. 
Qualitatively, the differences to the nominal case are very small 
concerning the structure of the mass-distance diagram. As for the 
KS tests, the largest difference is obtained in the mass distribu- 
tion (Sm — 47% instead of 95.6%). At first glance, this seems to 
be a large change and shows that the KS tests are quite sensitive 
and large differences can occur while "by eye" changes appear 
quite small. In many other cases however, the KS significance 
was found drop much more e.g. to < 0.1%. The other KS results 
are 53, 57 and 37% for S 3 -m, S d and S[Fe/H], respectively. 

5.9.2. Varying 

The current knowledge about type I migration which embedded 
low mass planets undergo is rather shaky. For example, Nelson 
& Papaloizou (2004J have shown that migration of low mass 
planets in a disc with magnetohydrodynamic turbulence occurs 
in the form of a random walk, and that averaged torques don't 
converge to well-defined values for low mass planets. Menou & 
Goodman (2003 ) have found a strong sensitivity of local migra- 
tion rates to the background disk structure like opacity transi- 
tions. Paardekooper & Mellema (2006) and Baruteau & Masset 
(2008) have also shown that depending on specific local ther- 
modynamic disk properties, the net torque can also be positive, 
i.e. causing outward migration. In addition, poorly ionized "dead 
zones" in the gaseous disk with a very low effective viscosity 
may also play an important role (Ida & Lin 2008b). For example, 
they allow much smaller planets to open a gap and to migrate in 
type II migration characterized by very low a values rather than 
in type I (Matsumura et al. 120071) . 

In summary, it seems increasingly likely that type I mi- 
gration rates found in laminar, isothermal disks (Tanaka et al. 
2002 ) might not be the whole story and that the true migration 
rate of small mass planets might be much more complicated. 
Nevertheless, short of a better description, we use the rates given 
by Tanaka's et al. (2002) and multiply them by an arbitrary ef- 
ficiency factor f\. We do not claim that this factor captures the 
true type I migration rate. It merely provides a convenient way 
to investigate the changes in the characteristics of the synthetic 
planet population with changing migration rate (see also Ida & 
Lin|2008a]l. 

We have synthesized three populations which differ from the 
nominal case with /j = 0.001 by the magnitude of the type I 
efficiency factor fi, for which we use values of 0.01, 0.1, and 
1.0. In Paper I, planetary formation tracks for two different fi 
cases were shown (0.001 and 0.1), and some effects on the the 



planetary populations were discussed. Figure. 10 shows the fi- 
nal positions of 20 000 planets in the a - M sin i plain for the 
population with /i=0.001, 0.01, 0.1 and 1.0. 

When generating initial conditions, we don't include plane- 
tary migration, as explained in Paper I. This implicitly means 
that we neglect the radial drift of the embryos at masses be- 
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Fig. 10. Semimajor axis a versus projected mass M sin i for four synthetic populations differing in the type I migration rate efficiency 
factor fi. In all cases, 20 000 synthetic planets are plotted. The dotted line is the feeding limit at a touc h. Planets migrating to this limit 
have been plotted at 0.1 AU. The fi = 0.001 case (top left), is the nominal case i.e. the plotted planets are simply a subset of those 
in fig. [2] panel A. 



low M em b,o = 0.6 M e . For low values of the efficiency factor 
fi (0.001-0.1), this is justified, as the type I migration rate in- 
creases linearly with planet mass (e.g. Alibert et al. 2005a]l, so 
that the extend of migration a seed would undergo while its mass 
is below M em b,o is on ly small. This can be seen in the formation 
tracks in Paper I. For fi = 1 .0 however this approximation is no 
longer correct, as illustrated by the results of Ida & Lin (2008a). 
In order to attenuate the artificial reduction of type I migration 
we have therefore synthesized the fi = 1 population with an ini- 
tial mass of the seed of 0.1 instead of 0.6 M m . 

At fi = 0.01 (top right panel), the synthetic population is 
in general similar as in the 0.001 case. A separation between 
the "main clump" and the "outer group" has however become 



visible, and the close-in very low mass planets (a < 1 AU, m < 4 
M e ), for which the effects of type I migration are the most severe 
(see also Ida & Lin 2008a) have started to migrate and grow 
beyond the isolation mass, and sometimes also have reached the 
feeding limit. 

At f\ — 0. 1 (bottom left panel) the effects of the type I mi- 
gration start to become clearly visible compared to the nominal 
case. The sub-population of "failed cores" is clearly reduced at 
large distances. Most embryos starting inside the iceline have 
fallen into the feeding limit, and only "failed cores" starting late 
and beyond the iceline are retained. Nevertheless, some features 
of the nominal populations like the "outer group" and the "hori- 
zontal branch" are visible also at this higher migration rate (by a 
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factor 100). The inner boundary of the "main clump" is shifted 
inwards and reaches to the inner boarder of the computational 
domain. A significant number of planets with masses up to 20 
Mq^_ have now migrated to the feeding limit. The presence of 
these massive planets is one of the reasons for which we re- 
jected high values of fi since these objects should have been 
easily detected. In fact, such planets have only been detected or- 
biting host star which are members of multiple stellar systems, 
but not around single host stars (Udry et al. 2003 ; Eggenberger 
et al. 120041 . 

The faster type I migration has also, as expected, conse- 
quences for the low mass planets reaching the feeding limit. 
In particular, with increasing value fi smaller and smaller mass 
planets migrate and eventually reach 0.1 AU, For example, 
namely 6 and 2 M e planets reach 0.1 AU for fi — 0.1 and 1.0, 
respectively. In the case of the nominal population only plan- 
ets with a mass larger than about 10 M e reach the inner bound- 
ary of the computational disk. In addition, the nominal model 
shows a depletion of planets in the region between the "failed 
cores" which virtually do not migrate, and those in the "hori- 
zontal branch" as discussed in §5.2 This depletion disappears at 
higher values of fi. 

Finally, we note that even for f\ - 1.0 giant planets can 
form (Thommes & Murray 2006), but only at a reduced number 
(tab. [5j. The sub-population of "failed cores" which was very 
numerous just beyond the iceline at small type I rates has al- 
most completely migrated inside 1-2 AU or even into the feeding 
limit. This has strongly populated a part of the a - M sin i plane 
(M sin i ~ 4 M B , a < 0.3 AU) which is depleted at low fi. 

The presence or absence of essentially gas-free, close-in icy 
or "ocean" planets (Leger et al. 2004) with masses below 10 
M m for which we can exclude that they have lost a significant 
primordial atmosphere could be used as a strong indicator of 
the efficiency of type I migration. Indeed we find in the model 
that planets which are (1) mainly icy (> 50 % of the accreted 
solids are icy), (2) which have a mass smaller than 10 M m , (3) 
a Me nv /M core that is much less than 0.1 and (4) which migrate 
from beyond the iceline to a < 0.1 AU only exist if f\ = 1. At 
fi = 0.001 we also find planets in the feeding limit which have 
a mainly icy core, but these planets have a mass of at least 20 
Me, and have more massive envelopes 0.1 < M env /M core < 1.0, 
i.e. they are of a Neptunian nature. Hence, the discovery in large 
numbers of low mass "ocean" planets would indicate that the 
type I migration as described in our model is not correct because 
the existence of such planets requires a high efficiency for this 
migration while the more massive and distant planets call for a 
sharply reduced rate. 

We have statistically compared all different fj populations 
with the observational comparison sample. Table [5] gives an 
overview of the results of the KS tests. 

The inspection of tab. [5] shows that both the results for 
fi = 0.001 and fi = 0.01 are all in at least fair agreement with 
the observations, with the nominal, slower migration case having 
somewhat better results for S a -M and S a , and the population with 
fi = 0.01 reproducing better the constraints derived for the "Hot" 
planets, and the metallicity distribution, S[F e /H]- At / : = 0.1 
the mass is still fairly well reproduced, but the significances for 
a - Msini and a are clearly reduced. At fi = 1.0 finally, all 
KS results have fallen to low or very low values, in particular 
for the semimajor axis distribution. The main reason for this is, 
as mentioned, that too many giant planets end up at distances < 
1 AU when the type I migration rate is high. Additionally, the 
period/semimajor axis gap between the "main clump" and the 
"outer group" which opens up at larger fi is not compatible with 



the observations either. The gap opens because the seeds that 
start inside the iceline in metal rich disks and which later become 
giant planets, now end up at smaller final semimajor axes than 
in the slow type I case, where they populate the region around 1 
AU. Therefore, the comparison of the synthetic populations with 
giant planets found by rather low precision (10 m/s) RV surveys 
constrain the effective type I migration rate to be < 1/100 of 
Tanaka's et al. J2002I I result, although also a reduction factor of 
< 1/10 can probably not be completely rejected. 

Table [5] shows that the detection probability P that does not 
include the planets that reach the feeding limit is nearly inde- 
pendent of fi between 0.001 and 0.1, and only decreases by a 
factor of about 3 for /i=l. It is then down at 3.1 %, which is a 
factor 2-3 lower than the observationally determined percentage. 
For the detection probability using the simpler FV05 bias (Pfv) 
the same behavior is found. This is due to the fact that higher 
migration rates have two partially compensating effects: While 
a high type I rate drains seeds into the feeding limit (which re- 
duces the detection probability), it also facilitates their growth 
because they get quicker into regions where new planetesimals 
can be accreted (which increases the detection probability). The 
results here show that for a large span of values for fi, the two 
effects approximately compensate each other. Only for the pop- 
ulation with fi = 1.0, the first effect dominates. Note that this 
prediction is in contrast with Ida & Lin (2008a) who have found 
that the fraction of stars with giant planets is a monotonically de- 
creasing function of the type I migration speed, unless the effect 
of an enhancement of the solid and gas surface density near the 
iceline due to a "dead zone" is included (Ida & Lin 2008b ). 

The detection probability where we have assumed that all 
embryos that reach the inner boarder of the computational disk 
become detectable, regardless of their mass as it would be 
the case for a "perfect" detection technique (Pwhot), increases 
strongly with fi, because the fraction of initial conditions which 
ultimately lead to a planet migrating to a t ouch also increases 
rapidly with fi: For the nominal population iVhot/Mnit is just 2.0 
%, but at fi = 1.0 45 % of all seed embryo fall into the feed- 
ing limit i.e. roughly 20 times more. The fraction of Hot Jupiters 
detectable with the FV05 bias i.e. massive planets in the feed- 
ing limit in contrast changes only by a factor ~ 5 with f\: It is 
0.4-0.6, 0.7-1.0,2.0-2.9, 1.9-3.3, for /^O.OOl, 0.01, 0.1 and 1.0, 
respectively, where the two values again assume a stopping dis- 
tance of 0.1 or 0.01 AU. This means that a high type I migration 
rate increases a lot the number of close-in low mass planets (Hot 
Super Earth and Hot Neptunes), but not that much the number of 
close-in giant (Hot Jupiters) planets, as we expect it. To quantify 
this, we study below the initial mass distribution of the "Hot" 
planets. 

We note finally that even if our results favor a low type I 
migration rate, a non-reduced type I migration rate does not pre- 
vent the formation of all giant planets, as mentioned before. This 
result could appear at odd with some of our former results (e.g. 
Alibert et al. I2005al) . where a reduced type I migration was a 
pre-requisite for the retention of giant planets. However, in the 
afore-mentioned calculations, the embryos were introduced in 
the protoplanetary disk at f start = 0, neglecting their formation 
time, which led to a phase of rapid migration during the early 
evolution phases of the disk. In the present models, such an ini- 
tial phase of type I migration is avoided, since embryos are intro- 
duced in the disk after a few Myr, when the disk surface density, 
and therefore the type I migration rate, have decreased. This is 
in line with the results of Thommes & Murray (2006). 
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Table 5. Statistical comparison of synthetic populations calculated with different type I migration efficiency factors f\. All other 
parameters are kept constant, except that for the f\ = 1.0 the seed mass is 0.1 instead of 0.6 M @ . The quantities are the same as in 
table [TJ and [3] Where possible, observational constraints are also given (Naef et al. 2005] FV05). 



Feature 


/: = 0.001 


0.01 


0.1 


1.0 


Obs. 




Detection probability w/o N ho , (P) [%] 


8.7 


9.4 


8.8 


3.1 


~5- 


10 


Detection probability w. iV hot (-P w hot) [%] 


10.7 


14.0 


28.7 


48.0 


> 7.3 


± 1.5 


Fraction of cases migrating to a t ouch (Miot/Winit) [%] 


2.0 


4.7 


19.9 


45.0 






Significance KS a - M sin i (S ^-m) [%] 


87.7 


35.2 


7.8 


0.3 






Significance KS M sin i (5m) [%] 


95.6 


93.6 


57.3 


5.0 






Significance KS a (S a ) [%] 


63.9 


18.4 


6.7 


0.2 






Significance KS [Fe/H] (5 [Fe / H ]) [%] 


21.7 


40.8 


28.0 


10.0 






Detect, prob. with FV05 bias excluding "Hot" planets (Pfv) [%] 


4.8 


5.3 


5.5 


2.2 


4.5 ± 


0.7 


Detect, prob. with FV05 bias including "Hot" planets (/Vv.whot) [%] 


5.2-5.4 


6.0-6.3 


7.5-8.4 


4.1-5.5 


> 5.5 


±0.8 


Fract. of "Hot" detectable planets of all FV05 detectable planets(FFv) [%] 


7.1-10.7 


11.4-16.4 


26.4-34.8 


47.1-60.5 


> 19 


t 6 



5.9.3. The IMF of the "Hot" planets 

From the mass-distance plots for the four type I migration rates, 
and the considerations made before, one already guesses that 
the mass histogram of the planets reaching the feeding limit at 
«0.1 AU must be significantly different for the four fi cases. As 
the sub-population of planets close to the parent star (including 
also very low mass objects) is the one that will become accessi- 
ble to observations (both high accuracy RV and transit surveys 
from space as CoRoT or Kepler) more quickly than planets fur- 
ther away, it is of interest to study the IMF of these planets. We 
again express the caveat that when comparing the synthetic mass 
spectra with the results of upcoming surveys, one must keep in 
mind the many uncertainties affecting the "Hot" planets men- 
tioned in Paper I, in particular those of small mass. For example 
we neglect subsequent evaporation, and only address one (in- 
ward migration) of several other possible formation channels for 
low mass, close-in planets like in-situ formation or shepherd- 
ing by giant planet migration. See Raymond et al. (2008) for an 
overview, and also Kennedy & Kenyon (2008]). 

Note also that here no observational bias has been applied 
and that the stellar mass is still fixed to 1 M Q . These points will 
be addressed in a later work. Nevertheless, the mass spectrum of 
the "Hot" planets could be a good tool to better understand type 
I migration. 

In fig. [TT]the distributions of M sin i of the "Hot" planets are 
plotted for f\ = 0.001 and 1 .0. The other two fi cases lie between 



these curves. The third curve in the figure is discussed in 5 5.9.5 



In both cases a peak at low masses is seen which corresponds 
to the Neptunian peak found in the overall (all semimajor axes) 
IMF (fig. [3j. It is the result of planets migrating inwards in the 
"horizontal branch". The exact location and the smallest mass 
inside the feeding limit is systematically varying with fi: The 
higher fi, the more the mass peak is shifted to lower masses and 
the more it broadens. The higher fi also, the lower the mass of 
the lowest mass planet in the "Hot" population. At f\ = 0.001, 
the peak of the distribution lies at about 20 M @ , and few planets 
smaller than ~ 10 M m have been brought in. For the fast type 
I migration case, the peak lies at about 5 M e , and the smallest 
planet brought into the feeding limit has a projected mass below 
1M». 



Figure 1 1 also shows that the ratio of massive Hot Jupiters to 



small Hot Neptune and Hot Super Earth planets changes system- 
atically with fi, as mentioned in the previous section. To study 
this, we separate the "Hot" planets in a high and a low mass bin, 
using as in table|2]a separating (projected) mass of 30 M e , as we 
have identified a local minimum of the overall IMF at approxi- 
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Fig. 11. Distributions of the projected mass of the "Hot" planets 
(i.e. those that migrate to the inner boarder of the computational 
disk), for three different populations: The solid line is the nom- 
inal population with fi = 0.001. The dotted line is for fi = 1.0. 
The dashed line finally is for fi-Q.QQl too, but with the effect 



of the rockline included (see §5.9.5 1. The mass spectra show the 
result of only one (inward migration) of several other possible 
formation scenarios for close-in planets like in-situ formation, 
which are not included in the model. 



mately this mass^j We then find that there are 2.0, 3.2, 5.1 and 
1 1 .4 times more planets in the low mass than in the high mass 
bin for /i=0.001, 0.01, 0.1 and 1.0. This clear dependence, as 
well as the dependence of the absolute number of close-in low 
mass planets on fi discussed in the previous section will soon 
be observationally determined. Recent observational results in- 
dicate that close-in low mass planets are very common (Mayor 
etal. 120091 



2 Also the IMF of the "Hot" planets has a local minimum between 
the planets which do respectively do not undergo gas runaway accretion 
and is thus bimodal. It lies for the nominal population at a slightly larger 
mass than in the overall IMF, namely at M ~ 40 M e . 
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5.9.4. Heavy element content of "Hot" planets 

The combined measurements of the mass and radius of planets 
transiting their host star allow, at least within certain limits, to 
determine the relative fraction of hydrogen and helium and of 
heavy elements in their interior, which can for example be used 
to deduce constraints on the formation of the transiting planet 
(Figueira et al. 12008b . 

As the number of transiting exoplanets is growing quickly, 
it is interesting to check if there are statistical correlations be- 
tween stellar properties and the planetary composition. Indeed, 
as shown by Guillot et al. (2006), and confirmed by Burrows et 
al. (J2007 ) and Guillot (2008), internal structure models indicate 
that there is a positive correlation between the total amount of 
heavy elements in a planet Mz tot and the host star metallicity, 
provided that some 'missing physics' (which can be a modifi- 
cation of the equation of state, a higher opacity or an additional 
energy source) are assumed to be at work in a similar way in all 
transiting giant planets. 

Our results for the composition of the synthetic planets of the 
nominal population (/J = 0.001) reaching the feeding limit are as 
follows: There is no correlation of [Fe/H] with the maximal to- 
tal mass (envelope and heavy elements) M or the maximal mass 
of accreted planetesimals Mh ea vy But, as shown by fig. 12 there 



150 



is a positive correlation of the total maximal mass of heavy ele- 
ments Mz,tot (mass of accreted planetesimals plus mass of heavy 
elements accreted with the gas). This is shown by an absence of 
empty circles representing "Hot" planets with M > 100M ffl (as a 
first order approximation of the detection bias towards massive 
planets) with a high M Ztot at low [Fe/H]. This is similar to the 
results obtained by Guillot (2008 ) shown with large black cir- 
cles for the scenario of an additional energy source deep in the 
planet's interior. These observed planets also have masses larger 
than 100 M m . A scenario with increased opacities gives qualita- 
tively similar results. The synthetic population also reproduces 
the general finding of internal structure modeling that some tran- 
siting planets contain high amounts (> 100M ffi ) of heavy ele- 
ments. Note that the host stars in the sample of Guillot (2008) 
have masses between 0.8 and 1.3 M Q , while we have a fixed 
M« = 1M . 

The models of Guillot (2008) assume a concentration of the 
heavy elements in the core surrounded by an envelope with so- 
lar composition. Our results of a higher Mz tot at high [Fe/H] is 
however due to higher amounts of heavy elements in the enve- 
lope, and not in the core. It could therefore seem questionable if 
a comparison can be made, as in general the repartition can have 
a important influence on the radius (Baraffe et al. [2008). The 
"Hot" synthetic giant planets of primary interest here have how- 
ever a Mz,iot/M < 0.3 (typically ~ 0.1 as Jupiter and Saturn) 
so that the influence of the repartition of the heavy elements 
only has minor influence on the radius predictions (Baraffe et 
al. 120081 . 



The small filled circles in fig. 12 represent mainly the abun- 
dant sub-population of Hot Neptunes discussed earlier. They 
typically have 0.7 < Mz,iot/M < 1. The only currently known 
transiting planet around a solar like star of this type (HAT-P-1 1 
b, Bakos et al. 2009, M Zitot ~ 23M e ) is also shown in the figure 
and falls well in this sub-population where observational biases 
are still very important. 

We have excluded in fig. 



12 the handful of extremely high 



mass planets (> lOM^) in the feeding limit which are in the 
nominal population quite well separated from the rest of the 
"Hot" giant planets (M < iM^). The reason is that due to their 
extreme mass, their a touc h is very large (several tenths of an AU) 
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Fig. 12. Total mass of heavy elements of "Hot" planets as a func- 
tion of [Fe/H]. Large empty circles are synthetic planets in the 
feeding limit with a total mass larger than 100 M©. Small filled 
circles are synthetic planets with a total mass below this limit. 
The large filled circles with error bars are taken from Guillot 
(2008) and show Mz >to t calculated for 16 currently known tran- 
siting exoplanets, assuming a kinetic energy mechanism. HAT- 
P-1 1-bwas added at [Fe/H]=0.31 and M Ztot ~ 23 M e (Bakos et 
al. 120091 . 



so that it is not clear if they can be classified as "Hot" planets. 
In any case, these planets which start their formation inside the 
iceline in extremely metal rich disks ([Fe/H]> 0.4 and Eo Z 800 
g/cm 2 ) are extreme in terms of composition also: They typically 
have a M heavy of 100-200 M e and M Z tot of 200-800 M e , reminis- 
cent of the internal composition deduced by Baraffe et al. (2008 ) 
for HD 147506 b (aka HAT-P-2b, Bakos et al. 120071 where a 
higher disk mass (due to the higher primary mass) could com- 
pensate the lower [Fe/H]. 



5.9.5. The "rockline effect" 



In figure 1 1 we have also plotted the mass histogram of the plan- 
ets in the feeding limit for another population which is identical 
to the nominal case (i.e. fi - 0.001) except that we have included 
the rockline (Paper I). The corresponding population is plotted 
in fig. [T3] It illustrates how populations synthesis can be used 
to study the global consequences of a given theoretical descrip- 
tion of some physical mechanism, and to reject it if it leads to a 
population that is in disagreement with the observed planets. 

In this simulation it is thus assumed that at distances where 
the disk midplane temperature is higher than 1600 K at t = 
(the moment where the disk evolution starts) no planetesimals 
exist. This means that inside the rockline a loc k, the solid surface 
density drops to zero. This affects both the generation of the ini- 
tial conditions, as well as the formation of some planets later on. 
For the initial conditions it means that the minimal distance from 
where embryos can start must be larger than the rockline. 

During the formation of the planets on the other hand, the 
following interesting mechanism is observed as predicted by 
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Fig. 13. Planetary population obtained assuming that no plan- 
etesimals exist inside the parts of the disk where the initial tem- 
perature is larger than 1600 K. The "rockline effect" leads to the 
formation of many Hot Jupiters, but almost completely elimi- 
nates Hot Neptunes, which leads to a completely different IMF 
of the "Hot" planets. 



Papaloizou & Terquem (1999 ): Subcritical planets of the "hori- 
zontal branch" migrating in planet dominated type II get to the 
rockline, where they suddenly become supercritical, as the solid 
accretion rate drops to zero (Pollack et al. 1996; Alibert et al. 
2005b). This leads to a rapid accretion of gas, and therefore to a 
completely different mass histogram of the "Hot" planets, with a 
strong peak at about 1-2 M^, but almost no Hot Neptunes. At the 
small distance of a few 0. 1 AU, the mass of the rapidly growing 
planet quickly overcomes the local disk mass (Paper I), so that 
the migration rate goes down into the increasingly slow, planet 
dominated mode, which prevents many planets from migrating 
further in. An inner hole in the solid surface density therefore 
leads automatically to a pile up of about Jupiter mass planets 
at a distance which lies somewhere inside the position of the 
rockline a 10C k as seen by the migrating planet. This is very well 
seen in fig. [13] where the "horizontal branch" bends upwards at 
a distance of about 0.4 AU, the location of the rockline for mean 
disk masses (Paper I). The "rockline effect" could therefore be 
a disk thermodynamics controlled, coupled formation and also 
stopping mechanism for the observed Hot Jupiter population. 
Note that the "rockline effect" leads to a mass-distance (anti- 
correlation (the mass decreases with increasing distance) qual- 
itatively reminiscent of the observed correlation pointed out by 
Mazeh et al. ( 120051 ). 

The difficulty is to specify the relevant location (for the plan- 
etary formation process) of the rockline, as the disk temperature 
profile so close to the star changes very rapidly at the beginning 
of the evolution of the disk. The rockline initially moves quickly 
inwards, at least in an a disk without irradiation like the one we 
employ in our model, so that defining the location of the rockline 
at t = 0, as it has been done for this test, likely overestimates the 
importance of the effect. Observations indeed indicate smaller 
inner dust disk truncation radii (Eisner et al. 2005), which would 



shift the "rockline effect" closer in. This overestimation is also 
supported by the fact that the population obtained in this way 
has so many more close-in Hot Jupiters than Hot Neptunes that 
it seems incompatible with the recent detections of several Hot 
Neptune and Super Earth planets (e.g. Udry et al. [2006). The 
HARPS high precision program indicates that there are many 
more low mass, close-in planets than Hot Jupiters (Mayor et al. 
2009 ), and also in the Keck data many such (not yet announced) 
candidates exist (Cumming et al. 2008). Indeed, when we again 
separate the "Hot" planets of this special population into two 
mass bins as in the previous section, we find that now there are 
15 times more massive planets than planets with a mass less than 
30 M m which is inconsistent with these observations. We there- 
fore conclude that modeling the possible absence of solids close 
to the star in the mentioned way is inappropriate, and must be 
replaced by a different theoretical description. 

In reality, planetesimals drift inwards due to gas drag and 
thus might follow the instantaneous location of the rockline, de- 
pending on the timescales of planetesimal drift versus rockline 
recession, an effect that we do not include in the model. The 
question whether the "rockline effect" can occur in nature de- 
pends thus on how the timescales of disk temperature evolution, 
planetesimal drag, and embryo arrival at small distances com- 
pare to each other. In any case, the mass histogram of the "Hot" 
planets (as the CoRoT and Kepler mission will provide) can tell 
us a lot about migration mechanisms, but also about the thermo- 
dynamic structure of the disk in proximity to the star. 

6. Predictions for extremely precise RV surveys 

In this second part of the result section, we return to the nom- 
inal synthetic population and use it to make predictions for the 
results of radial velocity surveys with an extreme precision. It 
also serves as an illustration how theoretical populations syn- 
thesis calculations can be used to estimate the impact of certain 
instrumental properties of a given detection technique. 

6.1. RV measurements at the 1 m/s and the 0. 1 level 

Up to this point, we have always used a radial velocity preci- 
sion equal 6rv,mc = 10 m/s and a survey length of trv.mc = 10 
years as an appropriate mean representation of past real extra- 
solar planet search programs. In the last few years however, de- 
tections mainly made with the HARPS spectrometer (Pepe et al. 
2004) made it clear that an RV precision of 1 m/s or even bet- 
ter has become possible, and that 0. 1 m/s should be possible too 
in future even though some issues as the intrinsic oscillations of 
stars must be overcome. 

These improvement have allowed to start exploring new, ex- 
citing parts of the planetary IMF in the Neptunian and Super 
Earth mass domain (e.g. Santos et al. 2004a; McArthur et al. 
120041 Lovis et al. 2006; Udry et al. [20071 Mayor et al. [2009 ). 
Here we study how the properties of the detectable extrasolar 
population change if 6rv goes down to 1 m/s or even 0.1 m/s. 

6.2. Detectable sub-population 

In fig. [14] we have plotted the a - M sin i diagram for a RV pre- 
cisions of the synthetic survey of 10 m/s (same as panel B in fig. 
[2k, 1 m/s and 0.1 m/s. The survey duration trv,mc is in all cases 
TO years. Note that the synthetic detection bias, which uses real 
data from the ELODIE survey, leads to a reduction of the detec- 
tion probability already at an induced radial velocity amplitude 
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that is larger than the detection threshold, but allows with a low 
probability also the detection of planets below it, which is par- 
ticularly well visible in the 0.1 m/s case, where many low mass 
planets with an amplitude of the order of 0.1 m/s exist. Fig. [I4| 
makes it obvious that the lower crv.mc, the closer the detectable 
population gets to the underlying full population. 

At each precision, new types of planets that become de- 
tectable put new constraints on planet formation theories. 

At 10 m/s, mainly Jovian planets are detected. Their fre- 
quency allows us to extract constraints on the timescales on 
which sufficiently massive cores must be built. Hot Jupiters and 
other giant planets at intermediate distances show that planetary 
migration is an important process which must be included in 
formation models. The "metallicity effect" among these Jovian 
planets is a strong indication that core accretion is the dominant 
giant gaseous planet production channel, as this trend is naturally 
reproduced by many formation models based on this paradigm 
( §5.7| i. The diversity of giant planets can be interpreted as the 
consequence of disk properties. Finally, the dependence of the 
overall maximal mass at a given semimajor axis constrains disk 
structures and stopping mechanisms, especially if a very large 
number of stars (Ge 2007) is observed over several years. 

At 1 m/s, the detection of many Neptunian planets in the up- 
per part of the "horizontal branch" out to several AU becomes 
possible, as well as the discovery of some Super Earth planets at 
distances below 1 AU. It also becomes possible to observation- 
ally determine the amount of depletion in the "planetary desert". 
This constrains the timescale of gas runaway accretion, and/or 
the exact efficiency with which gas accretion continues after a 
gap has been opened by a planet. The upper boundary of the 
"horizontal branch" indicates the mass where planets go into gas 
runaway accretion (30-40 M @ ). 

At 0. 1 m/s, almost the complete planetary population, espe- 
cially also Super Earth and terrestrial planets of > 1 M m at < 1 
AU will become detectable. The maybe most fundamental im- 
plication of the calculations presented here, namely that almost 
all stars with no (giant) planet detectable today should harbor 
low mass planets instead (tab. [2]) will be observationally put to 
the test. Additionally, at such a precision, the lower part of the 
population of subcritical Neptunian planets migrating inwards in 
the "horizontal branch" also becomes easily detectable, and will 
give interesting hints on the criteria for planets to migrate in disk 
dominated type II migration (planetary mass, disk scale height, 
gap opening criteria). It also becomes possible to determine if 
the planetary IMF is indeed trimodal or if this changed by late 
time growth. Finally, the abundance of planets with a mass of 
a few M @ inside 0.3 AU will help to constrain type I migration 
rates ( flggg) . 



Figure 15 shows the distribution of the projected masses 
of the full planetary population (corresponding to the IMF in 
fig. [3]) together with the mass histogram of the detectable sub- 
populations in fig. 14 at 10, 1 and 0.1 m/s. The distribution at 



10 m/s was compared to the observational comparison sample in 
fig.|4] At 10 m/s, the mass distribution has only one maximum, 
and a relatively simple shape. At 1 m/s in contrast, the minimum 
in the IMF at 30-40 M e and the maximum at 10-20 M e becomes 
visible, making the distribution bimodal. Interestingly, very re- 
cent discoveries at the precision of < 1 m/s (Mayor at al. 2009; 
Mayor & Udry 120081) also indicate such a shape of the mass 
function at the lowest masses we can currently detect. The com- 
parison with the results presented here must however be done 
very carefully, as effects of the primary mass (many very low 
mass planets orbit around M stars) and of planetary multiplic- 
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Fig. 15. Distribution of M sin i. Solid line: Full population. 
Dashed, dotted, dash-dotted lines: Detectable sub-populations of 
a 10 year RV survey at 10, 1 and 0. 1 m/s, respectively. All curves 
are normalized by the number of planets N synt in the full popula- 
tion. Note how the distribution becomes bimodal at a precision 
of 1 m/s. 



ity (most such very low mass planets are in multiple planetary 
systems) could be at work, too. 

Compared to the distribution at 10 m/s we see that the 1 
m/s distribution carries much more information (and thus con- 
straints) on the planetary formation process, like here on the ef- 

The 0. 1 m/s curve finally 



feet of gas runaway accretion {{ 5.3.1 



must be interpreted keeping in mind all the caveats concerning 
the incompleteness of the model at very low masses. This ap- 
plies for example to the predicted trimodal shape of the distribu- 
tion. A solid prediction of the model remains however that many 
Super Earth planets should be detected: At 0.1 m/s, planets with 
a mass less than 7 M ffl make up ~ 25% of all detected planets, 
and ~ 52% have a mass less than 30 Mm. 



Note that in fig. 15 the histogram for e.g. 0.1 m/s still lies 



below the distribution of the full population even at large masses. 
This is simply due to planets with a period longer than 10 years. 

6.3. Survey detection threshold 

After varying fRv.MC only, we can also ask how the overall detec- 
tion probability of the synthetic survey changes as a function of 



its duration trv.mc- In figure 16 the overall detection probabil 



ities without and with "Hot" planets (P and P W hot> respectively) 
are plotted as function of the synthetic survey's duration, again 
for 6rv,mc =10, 1 and 0.1 m/s. These fractions can be used as 
an estimation for the detection yield of real surveys, although 
that the one-embryo-per-disk limitation (54.1 1 is likely to turn 



the number of detectable planets artificially down. If we can use 
the Solar System with its four terrestrial planets as a guideline, 
then the effect should be particularly pronounced at low masses. 

In general, as expected, the longer the survey, and the bet- 
ter the instrumental precision, the higher the number of detected 



C. Mordasini et al.: Extrasolar planet population synthesis II 



23 



10 4 

1000 k- 



•g 100 

• I — I 

2 



10 

1 



?T0 'm/d 



:••-./^^*•l^tt^/^HW.■ff/";\^^*HKKBS , ■■■ — 



_i i i 




_l I I 



e^O.l'ni/V 




0.1 



0.1 



1 

a [AU] 




Fig. 14. The detectable sub-population for an assumed instrumental precision 6rv,mc of 10, 1 and 0.1 m/s. In all cases, the survey 
duration trv.mc is 10 years. It is likely that the population of Earth mass planets appearing at 0. 1 m/s is in reality more abundant than 
shown in the plot, where the number of such planets at small distances is just given by the fraction of disks with an isolation mass 
larger than M em b,o inside the iceline. Even so, the plot illustrates very well how many extrasolar planet remain to be discovered. 
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Fig. 16. Overall detection probabilities, i.e. the fraction of em- 
bryos that become detectable planets, for the synthetic survey 
as a function of its time duration t RVi mc f° r three different pre- 
cisions erv.mc- The solid line assumes that all planets reaching 
fltouch fall into the star (P), while the dotted line assumes that the 
get detectable, regardless of their actual mass at a t0U ch (^whot)- 
The curves for 0.1 m/s and to a lesser extent also the one for 1 
m/s must be considered as lower limits. 



planets. First, the detection probability increases quickly with 
ongoing observations, as the semimajor axis up to which planets 
are found increases. After some time (~10, 20 yr for £rv,mc=10, 
0.1 m/s, respectively) only few new detections occur and the 
curve flattens out. This can be understood by inspection of the 
a - Msini diagram: For eRv,MC = 10 m/s, it is mainly the fact 
that beyond ~ 3 AU, the maximal planet mass starts to decrease 



again. For e RV ,MC =0.1 m/s, the reason is the paucity of planets 
outside ~ 7 AU. 

The level at which the detection threshold stalls is ~ 9 - 1 1 % 
for 10 m/s. At 1 m/s precision, the detection rate flattens out at 
around 19 - 20 %. For the 0.1 m/s accuracy planets are found 
around 34 - 36 % of the observed stars after 20 years of ob- 
servation, and the saturation is not yet completely reached. This 
latter values, and to a smaller extent also the 1 m/s predictions, 
are lower limits only. The fact that we nevertheless find that it 
should be possible to detect at 0. 1 m/s extrasolar planets around 
at least every third star is promising for future extreme RV pre- 
cision instruments like ESPRESSO or CODEX (Pepe & Lovis 
2008; Cristiani et al. 120071 1. Note that the yield of such an ex- 
treme precision RV survey should already start at high values 
(~ 20% after 3 years). 

It is interesting to compare the synthetic results with some 
predictions based on observed data and extrapolation to longer 
periods: Marcy et al. (|2005) have extrapolated the total rate of 
occurrence of giant planets to be roughly 12% which is consis- 
tent with the number of long term radial velocity drifts in their 
data. Naef et al. (2005), by inverting the detection probability 
map of the ELODIE survey deduce a fraction of stars with plan- 
ets larger than 0.47 Mq^_ and a period smaller than 3900 d of 
7.3 + 1.5%. We find a somewhat larger fraction of synthetic 
planets fulfilling these criteria of 11.4 %. In a similar study, 
Cumming et al. (12008 ), find a occurrence rate of planets with 
masses 0.3 < Msim'/M^ < 15 of ~ 14 ± 2% inside 10 AU, 
using a flat extrapolation to larger semimajor axes. In the syn- 
thetic population, 9172/70 000=13. 1 % of the planets fall in this 
category, in good agreement. Outside 10 AU however, we pre- 
dict that the occurrence rate does no more raise a lot, as can 
deduced from the unbiased IMF (tab. |2]>: The three highest mass 
bins (M > 100 M ffi ) make up 14.3 % of all planets. The reason 
for this is that in our nominal model (a - M plot in fig. 13}, the 
highest mass planets can reach at ~ 10 AU is M max ~ 10 M^, 
but rapidly decreasing to ~ 100 M e at ~ 20 AU. In Paper I a 
number of possibilities were discussed that could influence the 
behavior of M max (a). Future observations, possibly with tech- 
niques especially apt in detecting massive planets at large dis- 
tances as astrometry (e.g. Unwin et al. 2008) or direct imaging 
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(e.g. SPHERE, Beuzit et al. 2007 ) will provide important in- 
sights and help to make our understanding of planet formation 
more complete. 

Developing reliable synthetic observational biases for these 
two techniques, as well as for transit measurements, which can 
then be coupled to synthetic planetary populations as demon- 
strated here for the RV technique will give very important stim- 
ulus to better understanding planet formation, as different de- 
tection techniques often have the ability to constrain different 
aspects of the planetary formation process in a complementary 
fashion. This is the case for the microlensing technique (Cassan 
& Kubas 2007 ) which has already demonstrated its ability to de- 
tect very low mass planets (Beaulieu et al. 2006 ) in a part of the 
mass-orbit plane which is not accessible to current RV surveys. 
Here we just note that these discoveries fit well into our pre- 
dictions that many very low mass planets exist at intermediate 
semimajor axes out to ~ 10 AU. 

7. Summary and conclusions 

We have used our extended core accretion model to synthe- 
size populations of extrasolar planets orbiting solar type stars. 
We have identified the subset of potentially detectable synthetic 
planets and compared them with a sub-sample of 32 actual ex- 
trasolar giant planets selected to satisfy the model assumptions. 

The subset of potentially detectable planets has been identi- 
fied using a RV detection bias model that takes into account the 
intrinsic instrumental precision as well as the duration of obser- 
vational surveys. To keep the model tractable, we have adopted 
a precision of 10 m/s and a survey length of 10 years as repre- 
sentative. 

We find that the synthetic survey has an overall detection 
probability of 8.7 %, in good agreement with the actual result. 
The total fraction of initial conditions that eventually lead to the 
formation of a planet more massive than 100 M @ is 14.3 %, again 
in good agreement with the observationally extrapolated total 
fraction of stars with a giant planet (Cumming et al. 2008 ). 

We have then made several Kolmogorov-Smirnov tests to 
compare the synthetic and observed distributions of some of the 
most important characteristics of the known planetary popula- 
tion. 

To compare whether observed and synthetic planets of sim- 
ilar masses are located at similar semimajor axes, we use a two 
dimensional KS test in the mass-orbit plane. We find a signifi- 
cance of 87.8 % that both the synthetic and the observed planets 
are drawn from the same parent distribution. In particular, the 
lack of massive planets at small semimajor axes in the observa- 
tions (Zucker & Mazeh 2002 ) is found as well in the models. 

The comparison of the mass distributions alone leads to a 
high KS significance of 95.6 %. The distribution has a peak at 
about 1-2 Mq^. The decrease at larger masses corresponds to a 
decrease in underlying, unbiased distribution, while the decrease 
towards smaller masses is a combined consequence of the de- 
crease of the underlying population and the detection bias. 

We have studied in detail this underlying planetary IMF, 
which is characterized by several minima and maxima with 
physical significance. The synthetic IMF has a global maximum 
at the lowest mass (~ 1 M @ ) that can occur in in our model, a 
next local maximum at about 15 M @ , and a third at about 1-2 
M<2j_. The first minimum at ~ 7 M 9 corresponds to the transi- 
tion between Super Earth and Neptunian planets. Note that since 
our PIMF corresponds to the mass distribution at the time the 
gaseous disk vanishes, it could be possible to fill-up this mini- 
mum by the merging of smaller mass planets. The next minimum 



at 30 M e is more robust and separates solid dominated planets 
from giant gaseous planets that undergo runaway gas accretion 
(Ida & Lin 2004a). It is very characteristic of the core accretion 
mechanism. At large masses (> 1000 M e ) the PIMF decrease 
with increasing mass. However, a long tail up to almost 40 
exists. 

The comparison of the semimajor axis leads to a KS signifi- 
cance of 63.9 %. We also find an raise of the distribution after a 
initially flat part in log(a), but in the model it occurs at about 2 
AU, rather than about 1 AU as in the observational data (Udry & 
Santos 2007). We interpret the fact that we reproduce the mass 
distribution better than the semimajor axis as an indication that 
we describe mass accretion in a less rudimentary way than mi- 
gration. It is also probable that the semimajor axis is more af- 
fected than the mass distribution by the consequences of N-body 
interactions after disk dispersal which are not included in the 
model (Thommes et al. 2008 ). 

At last, we have compared the metallicity distributions. We 
find that, similarly to observations (Santos et al. 2004a), the de- 
tectable synthetic planets are characterized by a metallicity dis- 
tribution that is shifted towards larger [Fe/H] by about 0. 1 dex 
compared to the initial distribution. However, a rather low sig- 
nificance of about 22 % is returned by the KS test. We suspect 
that a reason for this low significance could be connected to the 
extremely simple planetesimal disk model used in these calcula- 
tions. Despite this, the "metallicity effect", which is the observed 
increase of the detection probability with stellar metallicity is 
also present in the synthetic population. Even if the increase of 
the detection probability in the model seems to follow a slightly 
different (weaker) dependence on [Fe/H], we satisfy the obser- 
vational constraint within its error bars (Fischer & Valenti 2005 ). 

The next observational constrain we studied is the fraction 
of stars with a Hot Jupiter in orbit. Planets inside 0. 1 AU have 
been excluded from the previous quantitative comparisons, as 
the knowledge of the disk properties and processes close to the 
star are very uncertain. Assuming now for simplicity that the 
mass of the planets remains constant once they have reached 
the inner limit of our computational disk (~ 0.1 AU), we find 
a rate of occurrence of Hot Jupiters of about 0.4-0.6 %, com- 
patible with observations (Fressin et al. 120071 ). For these "Hot" 
synthetic planets, there is a positive correlation between the stel- 
lar metallicity and the maximal total amount of heavy elements 
they contain, in agreement with internal structure modeling for 
transiting exoplanets (Guillot et al. 120081 ). We have thus satis- 
fied constraints coming from transit surveys, too. Our popula- 
tion synthesis also shows that the formation timescales of giant 
planets (mean formation timescale of 3.5 Myr) are in agreement 
with observed disk lifetimes which is another important con- 
straint coming from the formation environment itself (Haisch et 
al. l200Tb . 

Our conclusions from these comparisons is that the core ac- 
cretion paradigm coupled with migration, even implemented in 
a much simplified manner, can generate a planet population with 
characteristics resembling the one of the actually detected ones. 
The crucial point for this conclusion is that the synthetic pop- 
ulation can reproduce a large number of different observational 
constraints simultaneously. 

This was not evident from the beginning as a large number 
of combinations of model parameters did not lead to acceptable 
results. For example, the efficiency of type I migration needed 
to be reduced considerably (10 to a 1000 times) from the linear 
rate by Tanaka et al. (2002) otherwise it would have been im- 
possible to reproduce the observed semimajor axis distribution 
(but see also Ida & Lin 2008b). An other example concerns the 
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possible absence of solids inside the rockline which would be 
an efficient mechanism to produce close-in giant planets, which 
could explain the observed pile-up. We believe that such results 
and constraints are the essence of the utility of the population 
synthesis efforts presented here. 

Finally, the nominal population allows us to make predic- 
tions about the planets which currently cannot be detected. As 
an illustrative example, we have shown the impact of improving 
the precision at which radial velocities can be measured. Our 
results indicate that the observed mass distribution becomes bi- 
modal at a RV precision of 1 m/s. The latest discoveries indeed 
point to such a feature (Mayor & Udry 2008 ). At a precision of 
even 0. 1 m/s, a large fraction of the underlying planetary pop- 
ulation will become detectable (planets will be found around at 
least 30-40 % of all FGK stars), and the observed mass function 
traces the characteristics of the underlying actual mass function. 
At a time when new high precision instruments are in the plan- 
ing, this could be used to quantitatively define the instrument 
requirements need to achieve a specific science goal. 
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